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ABSTRACT 


This thesis derives design tools for determining and improving performances of 
communication links that use M-QAM coherent demodulators associated with Viterbi 
soft-decision decoding (SDD) in Additive White Gaussian Noise (AWGN) and 
Nakagami-m channels. Performance analyses for 16-QAM, 64-QAM, QPSK and BPSK 
associated with up to three convolutional codes, including the one used by the IEEE 
802.11a standard and the duakk code, are presented as practical applications. The main 
tools relate to the analytical derivation of upper bounds of the probability of bit error (Pp) 
for any M-ary coherent demodulator followed by SDD, a methodology for improving an 
upper bound of P; tightening it to realistic data, and the obtaining of the 
specific B (d) spectrum for any convolutional code intended to operate with a certain M- 
symbol modulation. All derivations involve statistical considerations over the AWGN 
and Nakagami-m channels, as well as in-depth analyses of modulator constellations. The 
tools and models developed can provide great optimization to bandwidth limited system 
designs that require high data rates, especially the wireless ones. Consequently, they have 
great application to many fields of digital communications, such as cellular telephony, 


wireless networking, satellite links, ship-to-shore and ship-to-ship communications. 
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EXECUTIVE SUMMARY 


The increasing demand for high data rates in communication systems is forcing 
manufacturers and engineers to design links that will work close to the limit of the 
channel capacity. In other words, links that should operate at the highest data rate 
possible within the expected level of reliability. Within the military context, one possible 
application for these links would be the usage of the high-frequency (HF) spectrum to 
obtain over horizon communications that could reach more than 10 Mbps without 
satellite intermediation. Ships and troops could keep an intranet link operating in HF. On 
the other hand, commercial applications like cellular telephony and wireless networking 
have enormous interest in broadening their data rates consuming the minimum bandwidth 


possible. 


In this scenario, the M-quadrature amplitude modulation (M-QAM), a nor binary 
technique, performs a great role. It demonstrates high efficiency in bandwidth limited 
systems, but demands high power in order to obtain the required reliability. Due to this 
later characteristic, M-QAM should be always associated with an efficient forward-error 
correcting (FEC) code to operate in an acceptable signal to-noise ratio (SNR). At this 
point, one class of FEC that demonstrates excellent performance is the convolutional 
codes, especially when the receiver employs the most efficient decoding algorithm 
designed specifically for them, the Viterbi soft-decision decoding (SDD). Certainly, the 
association of M-QAM with SDD will be very important to meet the highest possible 


data rate requirements. 


The presence of M-QAM is more noticeable in wire communications than in 
wireless due to the nonexistence of the propagation obstacles, and the possibility of 
placing regenerators in the midway. Those propagation obstacles, such as multipath, 
shadowing and Doppler effect, end up producing fading in the receiver process, which 
decays considerably the SNR for wireless applications. However, the “big boom” of the 
cellular telephony and wireless networking are changing this scenario by forcing high 
data rates and, consequently, pulling M-QAM into the wireless world. As an example, 


XV 


most standards for the third-generation cellular system (3G) are being experimentally 


deployed based on M-QAM nowadays. 


Probably, due to being relative new into the wireless environment, the literature 
does not provide analytical models to determine performances of M-QAM with SDD up 
to this point. All analytical analyses involving SDD are done for binary modulation 


techniques, such as binary phase shift keying (BPSK). 


This thesis provides a deeper study in the major concepts involving M-QAM with 
SDD and derives tools for determining and improving performances of such systems. It 
considers the classical process of upper-bounding the probability of bit error for BPSK 
with SDD systems, and extends it, by means of a vector approach, to a generic model of 
norbinary coherent demodulators followed by SDD. Consequently, the model for M- 
QAM becones a trivial special case of the previous generic model. Regarding the 
channel corruption, the models are obtained for the additive white Gaussian noise 
(AWGN) scenario in a first consideration. As a second phase, in order to account for 
fading, the generic and M-QAM models are updated to consider the effects of frequency- 


nonselective small fading Nakagami-m channels. 


The models’ derivation resides on considering the worst case sequence of 
symbols being transmitted. However, the probability that all wrong decoded symbols in a 
sequence form a worst case is very small for large M’s. The thesis takes advantage of this 
fact and develops a technique to realistically improve performances for large M’s. As a 
result of such improvements, many upper bounds can be tightened to realistic views of 


the coding potential. 


Another interesting aspect arises when trying to utilize a convolutional code that 
is considered good for BPSK, into a M-QAM channel. Does it obtain good performance 
as well? Can a code perform well for every M-QAM, regardless of M? In order to answer 
those questions, a more intense attention is devoted to the B(d) spectrum and the code 
free distance, which are the most important convolutional code parameters responsible 
for determining the code performance. As a result, a numerical method for 


adapting B (d) spectrum to be used in a certain M-QAM channel is developed. Actually, 


XVI 


the method generates different B (d) spectra depending on M and based on the code’s 


transfer function matrix. 


As a practical application for all tools and models previously discussed, the thesis 
analyzes and plots the performances of four modulations: 16-QAM, 64-QAM, BPSK and 
QPSK operating in both AWGN and Nakagami fading channels. Each of the modulations 
is associated with the (2, 1) convolutional code used by IEEE 802.1 1a standard and the 
(2k, k) dual-k code. For the 16-QAM case, an extra code is also used due to 


demonstrating excellent performance with relative low data rate loss. 


Finally, the thesis summarizes some interesting results such as the range of good 
performances for the IEEE 802.1 1a code as well as its different B (d) spectra, validity of 
the improvement technique, formulas for upper bounds of M-QAM with SDD operating 
in either AWGN or Nakagami fading channels, and relative differences in performances 


between fading and non fading environments. 


The tools here developed can provide great optimization on designs of any M- 
QAM link and increase the efficiency of bandwidth- limited systems. Consequently, they 
could be considered extremely worthy for obtaining high data rates in many fields of 
digital communications, especially the wireless ones, such as cellular telephony, wireless 


networking, satellite links, ship-to-shore and ship-to-ship communications, etc. 
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I. INTRODUCTION 


A. BACKGROUND 

As the internet era imposes itself as being an absolute necessity in the modern 
world, it increases the demand on high data rate exchange for all existent digital 
communication systems in order to allow all colors, pictures, audio and video information 


to be transmitted. 


In the wireless world, many systems and standards, originally designed to carry 
just voice channels, which demand a meager 3.3 kHz of bandwidth, are being modified to 
meet the high data rate requirements. New standards are just being created and 
experimentally deployed, such as the third-generation cellular system (3G). This trend 
also reaches the military environment, which seeks the many options to provide high data 
links to mobile units such as ships, aircrafts and troops. Intranet services could be 
provided to an operating unit in a long distance communication by exploiting traditional 


channels formerly used to transmit just voice, such the high-frequency (HF) links. 


In the world of wire communication, the trend for higher data rates is much easier 
to resolve due to the nonexistence of fading, i.e., propagation obstacles that 
tremendously affect the wireless performance, such as multipath, shadowing and Doppler 
effect. As an example, a simple wire telephone line can be interfaced by a very affordable 
dia-up V.90 modem and be able to carry around 50 kbps with no changes to the carrier 
company systems [1]. If more bandwidth is needed, a very-high-digital subscriber line 
(VDSL) can provide a pair of wires with up to 51 Mbps downstream. Offices around the 
world are creating a shared information work environment by simply interconnecting 
their computers with local area network (LAN) cables and obtaining 100 Mbps of data 


rate. 


In this meantime, wireless communications was greatly expanded during the last 
decade due to great improvements in quality and affordability of the second- generation 
cellular systems (2G). The code division multiple-access (CDMA) technique was 


developed for 2G and is now regulated by the IS-95 standard, opening a new paradigm on 


sharing the electromagnetic medium among various users. Nothing would be more 
natural than to expect wireless to provide most of the useful services yet available by 
wire. For that purpose, the communication engineer community is devoting considerable 
efforts to overcome the propagation matters, especially within urban areas. The 
establishment of the IEEE 802.1 1a standard is a solid result of those efforts. It regulates 
wireless networking and provides data rates up to 54 Mbps when employing 64-QAM 


modulation and the orthogonal frequency-division multiplexing (OFDM) technique [2]. 


The new wireless standards previously cited are based on certain principles for 
maximizing the data rate in restricted bandwidth. One of them concerns the importance of 
utilizing a nombinary modulation technique, called M-QAM, to increase the channel 
capacity. Another concept highlights employing a forward error correction (FEC) coding 
scheme to compensate for the increase of the probability of bit error due to channel 


fading and M-QAM power inefficiency. 


The use of convolutional codes as FEC is known to provide excellent results. 
Additionally, their optimizations occur in the presence of Viterbi soft-decision decoding 
(SDD) [3]. However, the literature does not provide analytical models to determine 
performances of non-binary M-QAM modulations associated with SDD up to this point. 
All analytical analyses involving SDD are done for binary modulation techniques, such 
as binary phase shift keying (BPSK), and at most extended the results to quadrature phase 
shift keying (QPSK), considered as two orthogonal BPSK modulations occurring at the 
same time. Models to analyze performances of M-QAM with SDD would be very useful 
in the process of providing high data rates for wireless communications, especially for 


systems subjected to severe fading conditions. 


B. OBJECTIVE 
The objective of this thesis is to provide design tools and analytical models to 
analyze performances of M-QAM modulations using SDD either in the absence or 


presence of fading conditions. 


The main developed tools relate to: 


e the analytical derivation of upper bounds for the probability of bit error 
(Pb) when modeling any M-ary coherent demodulator followed by SDD in 
pure AWGN and in Nakagami fading channels, 


e the application of the previous models to obtain formulas for the special 
case of M-QAM modulations, 

e the obtainment of the Bid) spectrum for any convolutional code, 

e the necessary adaptation on the Bid) spectrum of a certain convolutional 
code originally designed for binary channels, to be used in nor-binary 
channels, 

e a technique that improves the performance analyses of M-QAM channels 


by realistically tightening the upper bounds of Pb, 


e the application of the previous tools for determining performance analyses 
of QPSK, 16-QAM and 64-QAM, as well as BPSK, when using dua-k 
and IEEE 802.1 1a convolutional codes in AWGN and Nakagami fading 
channels, and 


e an interesting option of convolutional code to associate with 16-QAM. 


C. THESIS OUTLINE 

Chapter I can be divided into two parts depending on the levels of information 
provided. The first part, compounded by Sections A, B and C, reviews the fundamentals 
on M-QAM modulation and convolutional codes. Most of the symbology utilized 
throughout the thesis is established there. The general formula for determining the 
probability of bit error associated with the use of FEC is explained together with the 
splitting of this process in two phases. The first looks for the B (d) spectrum that is a 
characteristic of the convolutional code, while the second searches for an upper bound for 
the probability of error in the SDD algorithm, P,(d), applied to a certain modulation 
constellation. The second part of Chapter II includes Section D and is devoted to 
explaining methods for determining the B (d) spectrum of any convolutional code. 
Accordingly, an analytical and a numerical methodology are shown together with a 


procedure for adapting B (d) spectra of convolutional codes originally designed to operate 


in binary channels, to be used in nor-binary channels. At the end of the chapter 
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the B (d) spectra of the code employed by the IEEE 802.11a standard are obtained for 
BPSK, QPSK, 16-QAM and 64-QAM as an exemplification of the exposed numerical 


methodology. The correspondent spectra are summarized in a table for quick access. 


In Chapter III, the upper bound of P, (d) for any M-ary coherent demodulator 
followed by a SDD is analytically derived for a pure AWGN scenario. This upper bound 
is further associated with the B (d) spectra previously obtained in Chapter II and applied 
in the special case of M-QAM rectangular constellations. The final result is a model for 
determining probabilities of bit error for any M-QAM modulation with SDD subjected 
only to AWGN. Furthermore, a technique that improves the upper bound of P, (d) trying 
to tighten it to realistic data is developed. As a practical application, the performance 
analyses of 16-QAM, 64-QAM, BPSK and QPSK, associated with up to three different 
convolutional codes are shown, in which the duak-k and the IEEE 802.11a code are 
included. The results of an excellent convolutional code to be used only with 16-QAM 
comes out of these analyses, which also take into account the technique of upper bound 


improvements discussed previously. 


Chapter IV discusses the basic concepts of Nakagami fading models and derives 
formulas for uncoded M-QAM in this particular environment. Following the same 
principles applied to the uncoded channels, the model for determining probabilities of bit 
error of nombinary modulations with SDD in AWGN obtained in Chapter III is expanded 
for Nakagami fading channels. The performance analyses of the four modulations in 
Chapter III are now repeated for the new model utilizing four fading figures m: 3, 2, 1 
and 0.5. The codes are the same as before and the technique of upper bound 


improvements is also applied. 


In Chapter V, a summary of the obtained design tools and models is presented 
together with the main results from their application to the specific cases of 16-QAM, 64- 
QAM, BPSK and QPSK operating with dua-k and the IEEE 802.11a convolutional 


codes. Also, it suggests the application of the results of this thesis in future work. 


The appendix finally provides an instance of how to apply the analytical method 


to derive a general formula for the B (d) spectra of duakk convolutional codes. 
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I. M-QAM WITH CONVOLUTIONAL CODES 


Convolutional codes are not easy to study unless all the basic concepts involving 
them, as well as the systems that surround them, are very well understood. This chapter 
recalls the fundamentals on M-QAM modulation, convolutional codes, and probability of 
bit errors. All these concepts will be extensively used throughout the thesis. The chapter 
also discusses in greater detail an important parameter of convolutional codes for 
determining performances: the B (d) spectrum. It demonstrates how to overcome some 


practical problems involving B (d) spectra, such as how to determine and adapt them to 


different modulations. 


A. M-QAM 

M-ary Quadrature Amplitude Modulation (M-QAM) is a nor-binary memoryless 
modulation technique in which one of M different symbols is transmitted per time using 
two orthogonal carriers (in quadrature). Each symbol is composed of g bits. The 


following formulas apply: 
= AB(D+ APCD, (2.1) 
q = log,M, (2.2) 


where 


a 2/T, 2nf.t), forO<ts< T, 
wo>| 5 cos(2T f. t) or 5 


0, otherwise 


0, otherwise, 


—./2/T, in(2 t) , forO<t<T 
vain={ 5 sin(27 f. t) or 


in which T, represents one symbol’s transmission time interval, f. is the carrier frequency, 


and A;, as well as Az, are the orthogonal carrier coefficients. 


The main advantage to using M-QAM is the gain in bit rate it provides. Notice 


that in the interval 7;, g bits are transmitted simultaneously and are embedded in only one 
) 


analog carrier. The final effect is an increasing of q times in the bit rate for M-QAM 
when compared with a binary modulation technique such as binary phase shift keying 
(BPSK). In other words, it increases the digital bit rate that can be transmitted within a 
limited bandwidth channel. One popular application of M-QAM is computer modems 
used for diakup connections when the user wants to maximize the communication bit rate 


over a restricted 3.3 kHz voice channel. 


This advantage does not come without a price, which is the increasing of the 
probability of bit error (P;) at a certain signal to-noise-ratio (SNR) in the channel when 
using M-QAM instead of a binary modulation. The larger the M, the worse the P» will be 
for the same SNR. In the end, the simplest way to keep P, at an adequate level required 
by the communication link is to add more power to the transmitted signal. The final result 
will be the enhancement of the SNR. Another factor to consider is the relative complexity 


of the transceivers that also grow proportionally with M. 


The M-QAM is completely defined by its constellation, which graphically 
represents the M possible symbols as dots on a Cartesian plot. The coordinates (x, y) of 
each constellation dot are the coefficients A; and Az of the correspondent symbol and are 
addressed by Equation (2.1). A popular constellation is one that is rectangular in shape, 
shown in Figure 2.1 for M equal to 16. This type is used by 2,400-bits/s V.22 bis 


computer modems [1]. 





Figure 2.1. |The 16-QAM Rectangular Constellation [After Ref. 1] 
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Even though all the concepts developed in this thesis are generic and could be 
applied to any M-QAM, a numerical focus will be given to four particular M's: 2, 4, 16 
and 64. These values have vast applications in numerous communication systems and 


define the BPSK, QPSK, 16-QAM and 64-QAM modulations, respectively. 


B. CONVOLUTIONAL CODES 

Another way to overcome the increasing of P, for a certain SNR when using M- 
QAM is to codify the information bits to be transmitted (channel coding) in such a way 
that only a small set of all possible received bit sequences are considered valid. The 
receipt of one invalid sequence indicates that an error has occurred in the channel. The 
decoder will then determine the closest valid sequence matching the received sequence 
and declare this valid sequence as the original sent instead of the one received. This 
process of encoding and decoding the information bits is called forward error correction 


(FEC) technique [1]. 


A special type of FEC commonly associated with M-QAM is one that combines 
previous input bits, stored in a memory, with the actual bits at the encoder input to 
generate the coded output bit sequence. Since new input bits constantly replace 
previously stored bits in a linear shift-register manner, the final result is a convolution 
operation between the input bit stream and a certain built-in function. This type of FEC is 


known as convolutional code. 


A convolutional code can be precisely described by some parameters. The main 
parameters are briefly depicted below but only in order to review basic concepts on the 
subject to be used throughout this thesis. Extensive detailed explanations are found in the 


reference literature. 


i Code Rate (R,) and Constraint Length (K) 


A (n, k) convolutional code causes the input bits to the encoder to be shifted into 
and along the shift register k bits at a time. The number of output bits for each k-bit input 


sequence is n bits [3]. Consequently, the code rate is defined by the formula 
Ro = kin. (2.3) 


Notice that n is always greater thank causing R; to be less than 1. Typical values for R. 
are between 1/4 and 7/8. Since R; is less than 1, it can be seen that the use of 
convolutional code reduces the final maximum bit rate capable of being transmitted in a 
bandpass communication channel. Also, notice that a small value for the code rate R, 
indicates a high degree of redundancy, which should provide more effective error control 


at the expense of increasing the bandwidth of the encoded signal [1]. 


The convolutional encoder memory consists of a shift register with K (k-bit) 
stages and n linear algebraic function generators (modulo-2 adders). The parameter K is 
called the constraint length of the convolution code [3]. Therefore, each one of the n 


output bits is affected by K sets of k input bits. 


2. Generator Polynomials and Transfer- Function Matrix 


A convenient way to represent the k-bit inputs and n-bit outputs along time is 
through polynomials in D, where the indeterminate D is interpreted as a delay operator. 
Each of the n modulo-2 adders is connected to the various shift register stages and is 
responsible for the generation of one of the n output bits. Consequently, it is possible to 
imagine a certain set of generator polynomials G\(D), organized in a matrix template, 
when multiplied by a polynomial representing the k-bit input, X(D), will lead to another 


polynomial representing the n-bit output, Y(D). The mathematical representation follows. 


Y(D) = X(D)G(D), (2.4) 


where 


Y(D) = [Y°@) YD) ... YD), 


X(D) = [ x) X(D) ... SD) |s 


and 
G°(D)  G(D)... ~=G"?(D) 
G°(D G(D ae GPO(WD 
Gp) - |& © 1D) aD) 
G&D) GD) se G27) 


The matrix G(D) is called the transfer-function matrix for the encoder [4]. The 
various generator polynomials G)(D) that compound G(D) are usually represented by 


their coefficients expressed in octal base. For example, (14)g means (001 100)2 and would 


represent Gi (D)=D°+D°. 


3. Trellis Diagram, Coded Symbols and Hamming Distance 


A very compact and convenient way to represent graphically all possible 
convolutional code outputs as a function of their inputs is by means of a trellis diagram. 
Each branch in the trellis expresses not only the n-bit output associated with the k-bit 
input, but also graphically indicates the correspondent encoder's internal state change. 
The set of successive branches defines a sequence of outputs also known as a path of 


code words. 


Figure 2.2 shows the n-bit output of the convolutional encoder feeding an M- 
QAM modulator that accepts qg bits per time. Considering n > q, those n-bits in N sets of 
q-bits can be grouped, where g is given by Equation (2.2) and N is essentially an integer 
number. On the other hand, if n < g, N sets of n bits will be necessary to make a group of 


q bits. In both cases, the qg bits are the ones that travel simultaneously through the channel 


in the form of a symbol, which is called a coded symbol. Therefore, each branch in the 
trellis can either represent N coded symbols (n > q) or only 1 coded symbol (” < q), 
depending on the situation. Furthermore, any path in the trellis, since it is composed of 


various branches, embraces many coded symbols. 


Coded Bit 
Sequence 


(n-bit stream) 







Information 
Symbol 


(k-bit stream) 


Conv. 
Encoder 





q-bit Shift 
Register 


(q bits) 


- Coded 
— M-QAM Symbol 


Figure 2.2. Block Diagram of a M-QAM Modulator Using a Convolutional Encoder 


The Hamming distance or simply distance between two paths y' and y” in the 
trellis, represented by d(yy"), is the number of coded symbols in which the two paths 
differ [4]. 


4. All-Zero Path, Minimum Free Distance (dic) and Transfer Function 


A convolutional encoder usually starts operating with its internal shift register 
completely cleared, 1.e., in the all-zero state. In this state, its response to an all- zero input 
will always be an all- zero output. In other words, the encoder will remain in the all- zero 
state forever until a nom-zero input occurs. The all-zero path is the horizontal line 
connecting all zero-state-dots (So) in the trellis diagram. It expresses the correct decision 
that a decoder should make when the transmission consists of a zero- information 
sequence. This is very useful in determining convolutional code performances. The 


method consists of constantly sending a zero- information sequence. In case a channel 
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error occurs, the probabilities and possibilities on the behavior of the decoder algorithm 


can be estimated by studying the paths that could leave and merge at the all- zero path. 


The minimum free distance, denoted dyj;ee, is the mintmum Hamming distance 
between all pairs of paths that start and stop at So [4]. The larger dye, the greater the 
possibility of correcting eventual channel errors. Usually the increase of dyrce is 
accompanied by an increase in the constraint length for the same code rate, which will 
result in enhancing the encoder complexity. One of the goals of achieving a good 
convolutional code is to obtain the largest dye for a same constraint length and same code 


rate. 


Another useful resource for studying convolutional codes that will be referenced 
later in this thesis is the transfer function. Basically, it shows a mathematical expression 
that reports to the contribution of all paths that necessarily leave and arrive at the state 
So. The expression below is an example [3] of a transfer function for a convolutional code 


withR. =1/3and K =3, 


N D° 


TDN = ———_——— 
( ) 1-2ND 


The meaning of the factors D and N can be better understood by expanding 7(D,N )in an 
infinite series as follows: 


N D® 


T(D, N) = ———— 
( ) 1-2ND 


= ND° + 2N*D* +4N°D” +...= Ya,N@?D", (2.5) 
d=6 


where 


_ 42"? (even d) 
¢ 0 (odd d). 


Looking at the series in Equation (2.5), the first term ND° informs that there is only one 
possible decoding path presenting one (exponent of NV) information symbol different than 


the zero- information symbol that corresponds to six (exponent of D) coded symbols 


different than the zero-coded-symbol. The second term in Equation (2.5), 2N *:D* , means 
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that there are two possible decoding paths presenting two (exponent of NV) information 
symbols different than the zero- information. symbol that correspond to eight (exponent of 
D) coded symbols different than the zero-coded-symbol. The interpretation of the other 


terms follows this same rule. 


The transfer function is obtained by building and solving a system of state 
equations based on the state diagram of the convolutional code circuit. Each possible 
state transition corresponds to one branch. For this reason, the definitions of D and N are 


related to branches. 


C. PROBABILITY OF BIT ERROR 

The most important figure of merit in a communication link is the probability of 
bit error for the information bits (Pp) as a function of the SNR. Every digital 
communication link has to comply with a certain maximum P, associated with a realistic 
SNR that could be present in the channel. The level of this threshold depends on how 


tolerant the receiving digital system is to errors . 


The formulas for determining P», versus SNR vary accordingly to the 
convolutional code, interference noise, demodulator and decoder algorithm employed. 
For convolutional codes, the optimum decoding involves a search through the trellis for 


the most probable sequence, what is achieved by means of the Viterbi algorithm [3]. 


Notice in Figure 2.2 that even though the convolutional encoder generates n bits 
at the output, the considered coded symbols are g bits long. Those g bits are those that 


really travel through the channel and are subjected to errors. 


A general formula for determining an upper bound on P; for M-QAM with a 


convolutional code can be derived as follows. Consider first: 
1 
c= yim : (2.6) 


where Py is the average probability of an information symbol (k-bit long), which feeds 


the convolutional encoder during the transmission, is retrieved in error at the end of the 
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reception process, and k is the number of information bits that compound an information 


symbol, i.e., the number of bits shifted into the encoder per time. 


The application of the concepts of union bound, first-event error and transfer 
function make it possible to study error probabilities for the Viterbi decoding of a 
convolutional code [3, 4]. One of the well known results of such studies is the general 


formula below that depicts an upper bound for Py . 


d= pee 


where B (d )is the total number of information symbols (k-bit long) different from the 


zero-symbol associated with all paths with distance d from the all- zero path, and P, (d) is 


the probability of choosing a wrong path in the trellis in a pairwise comparison of the all- 
zero path (correct path) with another path else (wrong path) with distance d from the all- 


zero path. 


Substituting Inequality (2.7) in Equation (2.6), the general formula for an upper 


bound of P, is obtained, 


P< 


: B(d)P,(d). (2.8) 


Notice by the definitions of Equation (2.7) that B (d) depends only on the 
convolutional code employed. It relates the number of information symbols in error with 
the number of coded symbols detected in error. Since B (d )is created by a convolution 
operation, its expression is mathematically feasible either by employing an analytical 
method or a numerical method as will be seen later in this chapter. 

On the other hand, P, (d) is completely independent on the convolutional code, 


being a result solely of the type of noise existent in the channel, the demodulation and the 
decoder algorithm. Its calculation is complex and usually done by upper bounding 


considerations over a pattern of an all- zero sequence transmission as explained earlier. 
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The independence between B (d )and P, (d ) makes it possible to split the problem 


of determining P» in two completely isolated phases that will be put together by Equation 


(2.8). The first phase consists of the analysis of the convoulutional code in order to 


determine B (d), whereas the second issue concerns the channel model as well as the 


receiving process in order to estimate a reliable upper bound for P, (d). 


D. DETERMINING THE PARAMETER B (d) 


There are two efficient methods for determining B (d). The first is analytically 


from the convolutional code transfer function, and the second is numerically from the 
generator polynomials. Both can become difficult for some codes, especially those that 


are good, requiring some simplification arguments or some knowledge on where to stop. 


if Obtaining B (d) from the Transfer Function 


This is an analytical process involving two steps [3]. The first is to expand the 


transfer functionT (D ,N ) of the convolutional code into an infinite series in powers of D 
and N. The second step consists of taking the partial derivative of this series with respect 
to N, and making N equal to | afterwards in order to obtain another series in powers of D 
only. The expression for B (d) becomes the equation that defines the coefficients of this 
last series as a function of d. 


The Appendix shows an application of this method for obtaining a useful formula 


for B (d) of the duatk convolutional code. 
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2. Obtaining B (d) from the Generator Polynomials 


The process of obtaining the transfer function for a convolutional code is not 
trivial. For large constraint lengths, the complexity of the process most of the times avoid 


reaching a final formula. Hopefully, by observing the abrupt tapering form of P, (d) as d 
increases, and that B (d) will be employed in Equation (2.8), it can usually be seen that 


only the first five terms of B (d ) are significant for determining P,,. Under these 
circumstances, a numerical method can be applied to the convolutional code core in order 


to precisely compute just the first few B (d)’s. 


The method uses the code generator polynomials for forcefully generating a large 
number of paths in the trellis beginning and ending at the zero state. These paths contain 
a relatively large number of branches, usually more than 3 times the constraint length, to 
cover a sufficient range of d’s. The process of path generating must follow a trend in an 
ascending order of d. The computation of the number of coded symbols different than the 
zero-symbol contained in each of those paths gives d. The number of information 


symbols different than the zero- information. symbol embedded in the information 
sequence related to the considered path provides the contribution for B (d )made by this 
specific path. If all those contributions for B (d) are continuously accumulated (summed) 
in a mult+dimensional vector, where each coordinate represents one d, after computing a 
large number of paths, the vector coordinates will express the B (d) spectrum. Since the 


path generation overall behavior follows an ascending order, it can be easily concluded 
that the accuracy of the first terms of B (d) increases as the number of generated paths 
becomes larger. Furthermore, those terms will freeze as soon as all paths that could 
contribute to their correspondent B (d) ’s have been generated. The process should stop 


when all terms of interest, usually the first three or five, are obtained. 
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The flowchart in Figure 2.3 shows the basic algorithm for computing the B (d ) as 


described above. It can be seen that the process first generates the informatior-bit 
sequence at block 2 and later the coded-bit sequence at block 5. Considering that the 
information sequence generation obeys an ascending order of numbers of information 
symbols in error, the correspondent coded sequences will present an overall behavior of 
ascending order of d. The overall word stands since among some successive paths, one 
path could have d less than its previous path due to the characteristics of convolutional 


codes. 


Notice the presence of a filter in block 4 in Figure 2.3 to avoid the codification of 
invalid information sequences. Recall that the Viterbi decoding process applied to an all- 
zero coded sequence is always trying to compare the metric of paths that could leave and 
merge with the all-zero path against the metric of the all- zero path itself. Thus, a path 
with a large number of branches that leaves and bounces off the all- zero path instead of 
merging with it should be rejected. In fact, this path, at the moment it bounces off the all- 
zero path, can be split into two distinct paths with smaller d, which was probably already 


computed earlier in the process. 
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Input of: 
° Transfer Function Matrix (G) 


° Parameter M of M-QAM 


Compute: 
Constraint length (K); g = log M 
Total length of information sequence 
including the zero-pad tail 
Vectorof B (d) ? O 


InfoSeq ? next information sequence 









CodeSeq ?  G (InfoSeq) 


d ?  #nomnzero coded sym. in CodeSeq 6 
InfoSymError ? # info symbol in error 


Accumulate InfoSymError in vector of B (d) 










Display results | 10 






Figure 2.3. Algorithm for Obtaining B (d) from the Generator Polynomials of a 
Convolutional Code 


ie 


Figure 2.4 illustrates an example of an allowed path with six branches. Notice that 


when it merges with the all-zero path; it remains there to the end. 


All- zero path 





Figure 2.4... Example of a Valid Path in the Trellis That Must Be Considered When 
Computing B (d) 


Conversely, Figure 2.5 shows a path that bounces off of the all-zero path and must 


be rejected. 


All- zero path Bounce point 





Figure 2.5. | Example of an Invalid Path in the Trellis That Must Be Rejected When 
Computing B (d) 
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Of course, the same computation principle leading to B (d)could be used to 


achieve other useful parameters such as ot (d), the total number of paths (code words) 
with distance d from the all-zero path. For this purpose, some extra computations should 


be added to block 7 in Figure 2.3. 


As can be seen, this forceful process of determining B (d )is based on a relatively 
large number of loops, dependent on the total length established for the information 
sequence. For instance, if someone works with a code that has k = 1 (binary input) and 
constraint length K =7, considering that the number of constraint lengths for determining 
information sequence length is fixed in 3, the algorithm would have to generate and 
check 2°°7 = 27! = 2,097,152 information sequences. It is easy to realize that this 


number grows exponentially with K and can lead to an expensive computation cost. 


Therefore, a large effort should be allocated to stop the process at block 8 in Figure 2.3. 


Some optimizations could be inserted in the process to avoid the generation of 
information sequences that a priori is known to yield to an invalid path. For example, if 
there were 12 branches instead of six in Figure 2.5, all combinations of paths beginning 
with their first six branches following the pattern indicated in the figure would be 
considered invalid. These combinations can be called a “family” and since the first 
member of this family, the path indicated in Figure 2.5, is considered discarded, the entire 
family can be a priori also discarded, and a considerable amount of computation 


resources can be saved. 


As a final observation, the minimum free distance (dj;e¢) of the code is also 
achievable at the end of the process. The B (d) spectrum (computed vector) will present 
zeros as the first terms, indicating that there are no paths with d = 1, 2,3,...,(d pe 1) 


that could leave and merge with the all-zero path. The first non-zero term in 


the B (d ) spectrum represents d,e- and can be easily obtained by searching in the final 


vector coordinates. This computation is performed in block 9 in Figure 2.3. 
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3. Adapting B (d) for Various Coded Symbol Lengths 


Sometimes the definition of B (d )is a little confusing and some details pass 
unnoticed. For example, as discussed previously, B (d)is a function of the convolutional 


code employed. However, the same convolutional code can generate different B (d) when 
attached to different M-QAM modulators. Notice that the distance d is related to the 
number of coded symbols different than the zero-symbol. As seen in Section II.B.4, the 


length of the coded symbol varies with M according to Equation (2.2). 


Consider a certain system of a convolutional encoder generating one fixed coded 
bit sequence as a response to a fixed information sequence at its input. If the option of 
connecting this system to one of the two M-QAM modulators presenting different M’s is 
available, the same coded bit sequence would be clustered in a different amount of gq bits 
in order to generate the coded symbols for each modulator considered. Consequently, 
different d’s would result, which would impact B (d ) directly, even though the 
information bit sequence and the convolutional encoder remains the same for the two 


systems. 


Very often a communication engineer knows that certain convolutional code 
works very well for a certain M-QAM channel and wants to improve the bandwidth 
utilization by employing a modulator with higher M and would like to estimate the 


performance impact that would result if the same convolutional code is used. However, 


the B (d) spectrum tabled for that code is specific for the original M-QAM channel. The 


question is how can a new B (d ) spectrum be achieved for the new M-QAM considered? 


Looking at Figure 2.3, block 6 is responsible for adapting information and coded 


sequences into information and coded symbols. It is this part of the algorithm that really 
takes into account the g bits of the M-QAM channel to reflect on the B (d ) estimation. 
Thus, the answer for the previously stated question is that if the generator polynomials 
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for the considered code is known, the engineer can simply use the algorithm depicted in 

Figure 2.3 provided with the generator polynomials coefficients and the new desired M. 
The flexibility on counting information and coded symbols for the same original 

sequences makes this process a very useful tool to map different B (d) spectra for the 


same convolutional code. Notice, however, that a convolutional code that works very 


well for a certain M-QAM channel may not perform when working with a different M. 


4. The B (d) Spectra for IEEE 802.11a Convolutional Code When 
Applied to BPSK, QPSK, 16-QAM and 64-QAM Modulators 


A program in Matlab was developed for this thesis based on the algorithm 
presented in Figure 2.3. One of the various codes checked by this program is the 
industry-standard convolutional code used by the IEEE 802.11a widely applied to 
wireless networks. The code has the following characteristics [5]: transfer function 
matrix G=[133 171],,R. =1/2, and K =7. Table 2.1 shows the B (d)’s spectra 
obtained for this code when used with four different M-QAM modulators. The parameter 


i denotes the increment that should be added on dyree to obtain d, i.e., d=d,,,, +i. 
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| iso | izt | is2 | i=3 | isa | 
Ps fe fa fo | 


ox feof |» fs |= | 
OAM insu | 19573 | 207085 


Table 2.1. The B (d)’s Spectra for the Industry-Standard Convolutional Code 
(R, =1/2) Used by the IEEE 802.11a When Applied to Different M-QAM Modulators. 





Notice in this case that the code rate is 1/2. Thus, the convolutional encoder can 
only be applied to M-QAM modulators presenting g equal to | or g as a multiple of 2, 
Le., {g =1 (BPSK or 2-QAM), 2 (QPSK or 4-QAM), 4 (16-QAM), 6 (64-QAM), 8 (256- 
QAM), ... }, through the arguments presented in Section II.B.4. Analogously, if the code 
rate were 1/3, allowed g should be {g =1 (BPSK or 2-QAM), 3 (8-QAM), 6 (64-QAM), 


9 (512-QAM), ... }. 


E. SUMMARY 

Chapter II reviewed some fundamentals on M-QAM modulation, convolutional 
codes, and probability of bit errors. The general formula for determining the probability 
of bit error associated with the use of FEC was highlighted, and special attention was 
given to issues related to the obtainment of B (d) spectrum. The most important part of 
Chapter II presented a numerical method for determining the B (d) spectrum of any 
convolutional code. As an application, the B (d) spectra of the IEEE 802.1 1a 
convolutional code were obtained for BPSK, QPSK, 16-QAM and 64-QAM. 
Additionally, a formula for determining B (d) spectra of dua-k codes was analytically 


derived. In the next chapter, all the concepts and results obtained thus far will be utilized 
22 


to obtain models, formulas and performance analyses of M-QAM with Viterbi soft- 


decision decoding in the AWGN scenario. 
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Il. PERFORMANCE ANALYSIS IN AN AWGN CHANNEL 


As explained in Section II.C, the most important figure of merit in a 
communication link is the probability of bit error for the information bits (Pp) as a 
function of the SNR. Therefore, the performance of every digital link can be estimated by 
plotting its P;, against the SNR. This chapter derives a model for determining an upper 
bound of P;, for any M-QAM modulation with Viterbi soft-decision decoding subjected 
only to AWGN. Furthermore, a technique that improves this upper bound will also be 
developed and applied to performance analyses involving 16-QAM, 64-QAM, BPSK and 
QPSK. 


The computation of P; is performed on the receiver part of the link, and its 
circuits are in an opposite order when compared with the transmitter as shown in Figure 


3.1. 








Analog 
Demodulated 

Signals 
Information 
Received Bit-sequence 

Signal Demodulator Decoder {0, 1} 
r(t) bee 
Figure 3.1. | Block Diagram of the Main Receiver Circuits. 


The parameter P; will depend on both the demodulator and the decoder. The 
dependence on the demodulator is related with how sensitive this circuit is to the channel 
noise, 1.e., how the noise associated with the signal at its input affects the demodulator 
outputs. On the other hand, the dependence on the decoder rests on the type of processing 
performed on its analog input (demodulator output) in order to decide which information 


bit sequence has been transmitted. 
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This thesis will concentrate on all the performance analyses considering the 
demodulation performed by a M-QAM coherent demodulator and the decoding done by a 
Viterbi soft-decision decoder (SDD). Notwithstanding this configuration requires more 
complexity from circuits and algorithms than when using a noncoherent demodulator 
and a hard-decision decoder, the results obtained are much better. The goal is to 
concentrate on the best performances since today’s technology provides relatively cheap 


resources to implement complex circuits on a large scale. 


A. VITERBI SOFT-DECISION DECODING 

The Viterbi soft-decision decoding (SDD) is a process of determining the 
information bit- sequence related to one or more sets of demodulated coded symbol 
outputs by using a maximum-likelihood sequence estimator (MLSE). The principle is to 
exploit the interdependence (memory) that exists among successive coded symbols, 
established by the convolutional encoder, in order to choose the most probable coded 
symbol that was actually transmitted given the reception of a certain signal corrupted by 


noise. 


Recalling that the trellis characterizes the memory in the transmitted signal, the 
MLSE algorithm searches among all possible paths through the trellis for a certain path 
that maximizes the probability that the received signals represent the coded symbols 
expressed by this path. For this purpose, the decoder associates a metric, called 
correlation metric, to each path. This metric is mathematically described for a certain 


path i, compounded by B branches, according to Ref. [3] as 
CMOS way log f.(RalC?), (3.1) 
j=l jal 


where f, (R, | Cc) represents the conditional probability density function (pdf) of the 


vector that expresses the received signals of the j-th branch (Rj) given that the coded 


symbols of the j-th branch of the i-th path ( Cc”) are considered to have been transmitted. 
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Notice by this definition that ig can be considered a branch metric causing the 
correlation metric to be the summation of all branch metrics in the path. 


After computing the correlation metric for all possible paths reaching a 
considered stage in the trellis (B branches ahead of some starting point), the MLSE 
simply chooses the path with the highest metric and discards all others. Considering that 
the number of paths to compute the metrics grows exponentially with the number of 
branches, some optimizations are done in order to save computation resources. For 
example, taking advantage of the fact that the correlation metric is cumulative, the 


algorithm needs only to keep track of a fixed number, 2(*~”” 


, of the highest metrics and 
their corresponding paths at a certain stage in order to compute the next stage metric. 


Those tracked paths are called survivors [3]. 


B. ADDITIVE WHITE GAUSSIAN NOISE CHANNEL 

The additive white Gaussian noise (AWGN) is a term that refers to the 
characteristics of the thermal noise, which is present in all communication systems and is 
the predominant noise source for many of them. It corrupts the signal in an additive 


fashion causing: 

r(t) = s(t) + n(t), (3.2) 
where r(t) is the total received signal, s(t) is the component due to the transmitted signal, 
and n(t) is the noise component. 


Thermal noise is considered white because its two-sided power spectral density 
G(f)=N,/2 is flat for all frequencies of interest [6]. It also presents a zero-mean 


Gaussian pdf depicted by 








fy(n) = —— diel 3.3) 


From this point forward the notation N (X, 6 ) will be used to denote a normal 
or Gaussian random variable (r.v.) X with mean X and variance o; . 
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The autocorrelation function of the AWGN is equal to 6 °5 (t) . Also, 


G,(f)=N,/2 is the Fourier transform of this autocorrelation, which causes: 
2 No 
o =—. (3.4) 
2 


A channel is designated as an AWGN channel when its impairments are limited to 
the degradation caused by the thermal noise. Particularly, the AWGN present in a signal 
is split within a coherent M-QAM demodulator into its two orthogonal channels. Such 
noise can be represented by a bidimensional vector of orthogonal AWGN components 


expressed by 
N, = [n,, no], Gs} 


where nx; and n,2 are independent- identical distributed (iid) Gaussian random variables 


with pdf’s given by Equation (3.3). 


c. DETERMINING THE PARAMETER P2(d) FOR M-QAM 

Recalling Section II.C, the type of noise in a channel affects only P,(d) . In this 
Section, a general formula for determining the upper bound of P, (d) for every M-ary 
modulation, whether orthogonal or nor-orthogonal, will be analytically derived for an 
AWGN channel and customized for the special case of M-QAM modulators that use a 
rectangular constellation. Furthermore, a way to improve this upper bound, 1.e., make it 
tighter, will be shown in an attempt to match the realistic communication characteristics 
for a particular link. The overall result will reflect directly on the predictions of the 


channel performance. 


iB Derivation of an Upper Bound of P2(d) for any M-ary Coherent 
Demodulator Followed by a Soft-Decision Viterbi Decoder 


In this Section the M-ary bidimensional case will be considered in a first analysis 
and then the results will be extended to a multi-dimensional modulation scheme. 
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Assume the block diagram in Figure 3.2. 


Y(t) 
$ q Soft- 
Decision (0, 1} 
r(t) bse (t) Viterbi 


Decoder 








Figure 3.2. Block Diagram of an M-QAM Coherent Demodulator Followed by a Soft- 
Decision Viterbi Decoder. 


The goal is to determine an upper bound for P, (d), which is the probability of 
choosing a wrong path in the trellis in a pairwise comparison of the all- zero path with 
another path that differs in d symbols from the all-zero path. For this purpose, the 


symbology that follows will be established. 


The vector C ,,, is the coded symbol related to the m'" symbol of the j‘" branch in 


the trellis. Recall from Section II.2.4 that the n-bit output of the convolutional code can 


be split in L sets of g-bit to adjust itself to the channel symbol length. This vector is 


and c, 


depicted by its two components as Ce, : Cina | , where c is 


jml can assume any 


discrete value mapped in the modulation constellation. 


The vector R jm tepresents the two outputs from the demodulator, 1.e., the Viterbi 
decoder inputs, corresponding to the m"™ symbol of the j" branch in the trellis. It is 


are shown in Figure 3.2. 


expressed as R ,,, = [ ee fas | , where r,,,, and r,,,, 


The vector N ,,,is the AWGN vector representation like the one in Equation (3.5) 


that affects the reception of the m'" symbol of the j branch in the trellis. Equation (3.2) 


can be written in a vector form as 
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Ry =C mp tN ju (3.6) 


Consider the following metric definition for the 7" branch of the i” path in the 


trellis: 
: = log SrA, Ri» R;, nse Ri }I {Cy Gy Cy . Cy }) 
= log fp{R,» m=1, 2,...N}I{C), m=1, 2,... N}) (3.7) 
= log freA Rin } | (Ce) 
where N is the number of coded symbols transmitted per information symbol. 


Since Kj is independent of Kjm for any m different than 1, it follows 


We AdRy tCCp. = TIARA TCC? (3.8) 


Considering that the reception of Rjz is only being affected by Cj, the following 


relationship stands 








PR VC Te Ss CR lo. (3.9) 
Substituting Equation (3.9) in (3.8), and Equation (3.8) in (3.7) yields 
ty Sloe] Lath, KC = Sloe x Gc, (3.10) 
m=l m=1 
However, due to the orthogonality of components | and 2 as well as the 
characteristics of the coherent detection employed, the conditional probability of 
Equation (3.10) can be derived as 
Fr Rim IC‘) = Fie Tina | er) FarCimn2 | eae 
= Cm = Cm)? = ima — Sina)” 
1 207 1 26° 
= ———— e — ‘ 
2710 * 2710 * 
which implies 
@ 1 = nee cin 
Fr(Rim 'Cin = ae (3.11) 


Substituting Equation (3.11) in (3.10) results in 





==il (iy 2 
N 1 =H Rj — CO 


a? = Y log 20 











1 (i) 2 
= WR, C2 ll 


iH} 
= 
— 
° 
gq 
nN 
+ 
— 
ie) 
gq 
— 
® 
—"” 
IM 


m=l 


Since N is a constant for the considered convolutional code, the above equation 


can be rewritten as follows 








N 
w=K,+K, PUR, -COu’, (3.12) 
m=l 
where 
i Re. sd 
K, = N log ne , and Noes paw ee 


As seen in Section III.A, the Soft-Decision Viterbi decoder computes an entire 
path metric for all paths that reach an intended stage in order to discard the lower metrics 
and keep the highest one. The correlation metric for this decoder is given by Equation 


(3.1) and rewritten below: 


B 
CM? =)". (3.13) 


Jal 


Substituting Equation (3.12) in (3.13) yields 


B N 
CM? = E(k IR,,-CO ‘ 


ae (3.14) 


BN . 
= BK,+K,Y\) UR, - Cin W? 
j=l m=l 
Considering the linear property of the convolutional codes, a sequence of all- zero 
symbols (Co ) can be assumed to be transmitted without any loss of generality. Following 


this transmission pattern, the Viterbi decoder should choose the all- zero symbol path, 


called path 0, as the correct path to be processed. 
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Suppose that, due to channel errors, the Viterbi decoder chooses an incorrect path 
1 instead of the correct path 0. Also assume that this wrong path 1 contains d erroneous 


symbols, i.e., d symbols different than Co. Therefore, an initial formula for P, (d) can be 


written as 
P(d)= Pr[ CM >CM na] = Pr[ CM —CM® > 0] (3.15) 


Substituting Equation (3.14) in (3.15) results in 


P,@) 


BON 
P| +KYY IR, —* 


j=l m=1 j=l m=l 


PES (UR,, —C2 I? -1R,,- CO I 20] 


j=l m= 


= 


In Equation (3.12) recall that Kz is a negative constant, which leads to 


j=l m= 


BN 
P,(d)= ap V2 (I Rin — CP IP=IR,, — CP 17) s 0. (3.16) 


Considering that a sequence of all-zero symbols is the one being transmitted, the 


following stands 


Ce =C,, Vjandm 
(3.17) 
R,, =C,+N yp: 


Replacing Equations (3.17) in (3.16) yields 


jm 


na=P SIN, +N,,-COW=-IG+N,,- Gi?) | 


j=l m= 
BN 
p> (1 (C® — Cy) — Ng l? -IN jp ll?) S | 


=r ES (ey -c > W-2(C%, ~C,).¥,) <0} (3.18) 


j=l m=1 


a2 


Recalling that C‘? = Ce =C,, for all combinations of j and m except for d 


jm 


symbols, Equation (3.18) can be rewritten as follows 


l=1 


ridy=P| ¥ (ic Clr = Gyn) so] (3.19) 


Expanding the noise vector in its orthogonal components in Equation (3.19) 


results in 


d 
P,(d)= Py (UC — C1? -2(€ = €4,) 4-22 - cama) S0| (3.20) 


[=1 


At this point, notice that P,(d) will be very dependent on the wrong symbol 


sequence CY , for/= 1, 2, ... d. However, for a given sequence of ce, Le., path 1, it can 
be seen that the left side of the inequality in Equation (3.20) is a Gaussian random 
variable due to nj and nj being zero-mean Gaussian r.v.'s and the other terms being 
deterministic values. This random variable will be called X, which leads to 
X=N (X, of )=N [S, Eo: } (3.21) 
= 


7 i=l 
where 

Xo Sher S71. il, (3.22) 
and 


57, =4( (cP -— cy) +(C2 - cy)? Jo? =411C — CI?” 


=4X0°. (3.23) 


(1)> 


Here an upper bound can be established by assuming that all C;"’’s, for every / 


are equal to a worst-case symbol, Cw, in such a way that P, (d) becomes the highest 


possible. 
If Cwis a constant vector for every /, then Equation (3.21) can be rewritten as 


X= NX» of) = NG; dos): (3.24) 
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The problem now is to find a symbol Cwdifferent than Co among all possibilities 
within the modulation constellation that could be the worst-case symbol. For this 


purpose, a useful theorem should be derived involving normal distributions. 


Theorem: consider two Gaussian r.v.'s X and Y 
N (X,o0%) and N (¥,o;), 
then 
Pr[ X <0]> Pr[¥ <0] eX, <Yo,. 
Proof: 


Pr[ X < 0]>Pr[Y <0] 








Since the Q-function is a descending monotonic function, in order for the above 


inequality to be true, the following stands 


Notice that the reverse also works, which concludes the theorem’s proof. 


By the above theorem, it can be stated that, within a set of Gaussian r.v.'s with 
different means {X pt=L25.2:,N } , the one that has the lowest mean will comply with the 


following property 
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Pr[ X, <0]>Pr[X,,, <0], where aD Oa Ee 


except if one of the other r.v.'s with a higher mean presents such a large standard 


deviation that makes 
O,> (3.25) 


which would hurt the theorem. 


Therefore, C wis thus chosen to make the mean in Equation (3.24) the lowest 
possible and then checked if there is any other remaining choice that would make the 


Inequality (3.25) true. 


Looking at Equation (3.22), it can be seen that the minimum mean in Equation 


(3.24) is obtained by choosing a Cw that makes 


IIC,, — C,I?=Dy 


O,min ? 


(3.26) 


where Do, min 1s the Euclidean distance between Co and its closest neighbors in the 
constellation. Notice that there is at least one Cw for every constellation, but in some 
cases, there are more than one. However, for those cases, any closest neighbor chosen to 


be Cw will yield the same statistical results (mean and variance) because of Equations 
(3.22) and (3.23). 


Substituting Equation (3.26) in Equations (3.22) and (3.23) yields 


Xy = Dios > (3.27) 
and 
Sxy =4D) min O (3.28) 


The next step is to determine if there is another possible choice of Cy, called Cww, 


which would comply with the following relationships 


Xp OX yp Md: Xi O ye, <> Xi Ove (3.29) 


and could make the exception pointed by inequality (3.25) true. 
With Equations (3.22) and (3.23), the second inequality at (3.29) becomes 
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Kyi kee GO? SX ald Xe O°. 
Since both X,, and X,, are positive values due to Equation (3.22), it is possible 
to square both sides of the above inequality and obtain 
46 °X2 X,,, > 40°X2, X, 
Rie OD grap 
This latter result hurts the hypothesis formulated by the first inequality in 


Equation (3.29). It means that C ww does not really exist and the choice of Cwas the 


closest symbol to Co can be definitely considered the one that will conduct to an upper 


bound for P, (d) in Equation (3.20). 


Finally, it should be recalled that the left side of the inequality in Equation (3.20) 
represents a Gaussian random variable defined by Equation (3.24). Thus, an upper bound 


for P,(d) can be set as 


P,(d) < Pr[X <0] 


_d¥ y 3.30 
B(d)< 1-0 0-dX,, La ax. (3.30) 


do yy do... 


Substituting Equation (3.27) and (3.28) in (3.30) yields 








d Dy a” Dia dD; 


0,min 


BG): <0) | 9 faa 
. ”| Tea Bi, o° 4d Dy min O 40° 


which implies 





Pd) < Q|,}-——™ |, (3.31) 


where the last step stands due to o* = N,,/2, expressed in Equation (3.4). 


Notice that Do, min can vary depending on the Co chosen to be the zero-symbol in 


the constellation. Therefore, in order to guarantee that Inequality (3.31) will provide the 
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real upper bound, Co should be chosen to be the one that yields Do, min = Dnin , where Dinin 


is the minimum Euclidean distance in the constellation. Thus, the final formula to 


P,(d)< | ee } (3.32) 


Due to the orthogonality achieved by coherent demodulators, all the assumptions 


determine P, (d) becomes 





thus far are also valid for multi-dimensional M-ary modulations such as M-FSK, as well 
as unidimensional ones such as BPSK. Consequently, P,(d) found in Inequality (3.32) 


can be considered a general upper bound to be applied to every M-ary coherent 


demodulator followed by a soft-decision Viterbi decoder. 
2, Upper Bound of P2(d) for M-QAM 


In this Section, Equation (3.32) will be customized for the cases of BPSK and 
QPSK, special cases of M-QAM in which M = 2 and M = 4, respectively, as well as two 
other M-QAM coherent demodulators, 16-QAM and 64-QAM. All of them will be 
considered to be using rectangular constellations. Notice that Equation (3.32) can be 
applied to every type of constellation. The particularization for the rectangular one results 


since this type is the most frequently used in practice [3]. 


For rectangular constellations, regardless of M, Dyin is always equal to a constant 


2A, where A is one unit of the grid represented in Figure 2.1. Hence, it stands 


2 
Dain 


=4A’. (3.33) 


Applying Equation (3.33) in (3.32) leads to 


d(4A’) | 2dA } 





B(4)s@) J—— Fel (3.34) 


2N, 
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Equation (3.34) will be used to derive the well-known formula for a BPSK 


demodulator followed by a soft-decision Viterbi decoder [3]: 


P,() suse $ O(f21,Red J. 


Assume the following symbology: 
e¢ %s is the average energy of one coded symbol €,, /N, (average SNR for one coded 
symbol), 
e 2?» 1s the average energy of one coded bit / No (average SNR for one coded bit), and 


e is the average energy of one information bit / No (average SNR for one 
information bit). 


As regards these parameters, it can be stated that: 


€ 
teres (3.35) 
No 
Ves = DY cb > (3.36) 
and 
Vo =ReYs- (3.37) 


Substituting Equations (3.37) in (3.36) and the result into Equation (3.35) yields 
€., = Re, No- (3.38) 


For BPSK, considering each coded symbol to be equiprobable to occur in the 


channel, the energy per coded symbol can be expressed by: 


eg Seley | -y Pr[symbol C,]e,, = > (2.47) =A’. (3.39) 


i=l 


Replacing the results of Equations (3.38) and (3.39) in (3.34) leads to 


2de 2d N 
ria) so| Hee Jol a [adgR, )- (3.40) 
0 0 
Recalling that g = 1 bit for BPSK, Equation (3.40) becomes 


P,(d) spsy < O(af2dRe¥, J. (3.41) 


38 


This latter result is exactly the one expected to be found. 


For the QPSK, the energy per coded symbol can be written as: 


e.,=Ele,, |= ¥  Pr[symbol CJe,, =4 (24 )=2A°. (3.42) 


i=l 


Replacing the results of Equations (3.38) and (3.42) in (3.34) leads to 






2d £s 


< -o ARNG) (aR) (3.43) 


0 





Recalling that g = 2 bits for QPSK, Inequality (3.43) assumes its final format as 


P,(d) spse <Q(xf2dRc%, J. (3.44) 


Notice that this result is identical to the one obtained for BPSK. 


Now is the time to see how the upper bound of P, (d) for the 16-QAM behaves. 
The constellation is shown in Figure 2.1. Notice that Equations (3.34) and (3.6) are the 


same, even though qg changes to g = log, 16 = 4 bits. 


Using the same arguments as for BPSK, the following computation stands 
16 1 
e,,=Ele,, |=) Pr[C,]e., = 404° +10A7+10A? +18A’)=10A?. (3.45) 
i=l 


Replacing Equation (3.45) in (3.34) causes 


(3.46) 





Finally, substituting Equation (3.38) in (3.46) yields 
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dqR YN, 
P, C3 eee <0 arash, | : 
0 


BA) cca <0) fate } (3.47) 


For the 64-QAM, the development is the same but g = log, 64 = 6bits, leading to 


which results in 


€,, =E|€,, ]-yPr [C,]e,, =424”. (3.48) 


i=] 


Replacing Equation (3.48) in (3.34) results in 


(3.49) 





Therefore, by substituting Equation (3.38) in (3.46), the final formula becomes 


dqR.-Y,,N 6dR-Y, 
R(da ow 0 | Ral Jeol at } 
[2 
P,(@)64-oam <0) ZAR, } (3.50) 


3. Improvements on the Upper Bound of P2(d) 


which yields 


The core of the analytical development of the latter upper bound for P, (d ) resides 
in Equation (3.21). The worst-case scenario for the result of this equation was assumed. 
However, some questions can be posed. Is the worst-case assumption too conservative for 
a certain application in terms of leading to a loose upper bound? Can different hypotheses 


other than the worst case be assumed and still have a reliable but tighter upper bound? 
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How far can this process go? The answers to these questions are the subject of this 


Section. 


Rewrite Equation (3.21) in a different fashion by replacing Equations (3.22) and 
(3.23) in it. 


d d 
fe@ = N [Enc ci soc —Gil} (3.51) 


[=I 1=1 


The worst case was shown to happen when all d symbols are the ones closest to 


Co in the constellation. In this case I|C; —C, ll *= D;,,, which led the summation in 


min ? 


Equation (3.51) to be the smallest possible, d D*,., and the final result to be Equation 


(32): 

Notice, however, that there are many symbol combinations for filling all d slots 
when d = 3 and M = 16. The probability of all symbols being the closest ones to Co tends 
to be very small for those cases. Even though small these special symbol combinations 
must be considered for the worst-case scenario, in which 100% of all possible 

d 
combinations result in y II (Ong =€, ll Pd dD... . Consider all combinations that lead to 
l=1 
the second smallest possible value for the latter summation instead of the smallest one 
(worst case). The percentage of the total d-symbol combinations this second smallest 


value will cover can be expressed as 


P.,, =1—Pr[all d symbols have distance D,,,, to C, | 


cov 


(3.52) 
=1-Pr[ IC)” - Gll’=D;,, for! =1,2,...,d]=1-P.,,. 


If the worst-case probability P,,- for a certain link is very small, the second 
smallest possible value for the summation could be assumed as the worst case without 


significant losses in generality. 


Consider the definition of an average distance, D,y, dependent on the specific 


combination of all d symbols chosen to be the wrong path, given by: 


4] 


d 
Ver aCe Sap; 


l=1 
2 we ‘ ay 2 
Diy = Dl COSC: (3.53) 
1=1 


Notice, in the worst case that D*,, = D:,,, and for all other cases, D*,, > D-.,. The 


idea of making a different choice for the worst-case sequence impacts directly in 
Equation (3.32) in which D-,, will be replaced by D%,,. Even if D%,, is just a little bit 


2 


greater than D-,,, due to the Q-function, the parameter P, (d) will be considerably 


reduced. 


The bigger the parameters d and M, the smaller the worst-case probability Py. in 
Equation (3.52) will be, and consequently, the closer to 100% the P.o, will become, and 


the closer to D? 


~., the parameter D+, will reach. Recall that d is at least equal to djree, 
which usually is greater than3 for good convolutional codes. Summarizing the issue, the 


main goal to improve P,(d) would be to determine an acceptable D%,, for d= djree that 


can replace D°.. in Equation (3.32). Even though the effort begins by considering 


n 


d = dyree , aS Will be seen later, the same Diy found for dj;ee works rather well for the 


other terms of P, (d) with d greater than dyjee. 


Applying Equation (3.53) for all possible d-symbol combinations, it can be seen 


that D%,, lies on a finite and discrete domain of values. Therefore, the best way to study 


2 


the possible choices of D%,, is to plot its discrete probability density function (dpdf) for a 
certain dyree. Notice, however, that the dpdf will be very dependent on the symbol in the 
constellation chosen to be Co. Figure 3.3 shows a 16-QAM constellation with two 
different possible choices of Co. The circled numbers indicate the square of the Euclidean 
distance, D’ , normalized by the grid unit A’ for every symbol different than the 
considered Cg. For every / in Equation (3.53), each correspondent term |IC\”) — C, Il? 
can assume any value of D’ represented in Figure 3.3. Notoriously, the numbers for 


configuration (b) are much greater than (a). 
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(a) - Co is one of the most internal (b) - Co is one of the corner symbols 
symbols 


Figure 3.3. Square of the Normalized Euclidean Distance (D*/A’) Related to the 16- 
QAM Rectangular Constellation Symbols for 2 Different Choices of Co. 


The choices of Co that provide larger values for IIC\” — C,|I* will lead to greater 
Diy and, consequently, will achieve better results for P,(d), i.e., a smaller P,(d). In 
order to be conservative, this study will focus on rectangular constellations choosing Co 
to be the one resulting in the smallest D‘,, and the largest P,(d), which is the worst 
choice available in terms of P,(d). For any rectangular constellation, it is easy to see that 


the conservative choice of Co falls on any one of the most internal symbols because they 
present the largest number of neighbors with Euclidean distance equal to Dmin. The option 


(a) in Figure 3.3 represents this choice. 


Consider a practical example. Suppose the analysis of D‘,, for a 16-QAM 
rectangular constellation is associated with a dual4 convolutional code, i.e., duakk in 
whichk = 4 to match the 4-bit coded symbols of the 16-QAM. In the Appendix it is 
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shown that djree = 4 for dua-k convolutional codes. Thus, by computing the dpdf of D‘%,, 


for all possible sequences of 4 symbols in which each symbol can be any of the 15 


symbols different than Co, the plot shown in Figure 3.4 can be obtained. The set of 


different values for Diy , 1.e., the horizontal axis, are again normalized by the factor A’. 
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Figure 3.4. Dpdf of Diy for a 16-QAM Rectangular Constellation, d = 4, and 
Conservative Choice of Co. 


Even though the number of possible combinations for a sequence of 4 symbols is 
equal to 15* = 50625 , due to Equation (3.53), those combinations result in only 26 


different values for D‘,, . The point where D‘,, / A’ =4 represents the probability of one 


worst-case sequence ( D‘,, = D-..) occurring. This value can also be determined 


min 


analytically by considering the existence of 4 different closest symbols to Co in the 
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constellation. Thus, there are 4* possible sequences that will result from D’, = D-,,. The 


min 


probability that one of these sequences occurs is given by 


4 


= ca = 0.00506 = 0.5%. 
15 


we 


Substituting this result in Equation (3.52), it can be stated that the choice of 
Diy =5A° instead of D*,, =4A° would be valid for 99.5% of all possible transmissions 
in the channel and could lead to d = dee = 4. What about the other values of P, (d) used 
for the calculation of P, in which d > djee? Since their d is different from dfee , should a 
specific Diy be determined for each d or can the the value obtained for dy-c¢ be applied to 


all d’s? 


Notice that when d increases, the domain of the possible D%,, also increases and 
all probabilities of the previous dpdf will be diluted along this new longer domain. 
Consequently, by using a certain threshold determined for dye in a path presenting 
d > dree, the probability of obtaining Diy greater than this threshold is greater than this 
equivalent probability when using the same threshold for a path in which d = dee. In 
other words, a threshold that is good for dye will be even better for d > dfec. To 
summarize, the analytical conclusion of the 16-QAM exampk can be restated as: “... the 


choice of D‘,, =5A” instead of D,, =4A* would be valid for at least 99.5% of all 


transmissions in the channel’’. 


Figure 3.5 shows the difference in performances for the discussed example. The 


curve for Peo, = 99.5% was determined following the equation 








P =7 y B(d)PA(4))o 1 5 (3.54) 
d=d free 
However, 
me d (54°) d(5A°)e,, 
P,(d) _ =O AV =Q =O 
Dhig= Hy =5 A 2N, 2N, 28 IN 6 


d(5A°) d(5A°) 
Doin= Dey=5 & = Q Ie Ves = ys GRY, 


cs cs 








P,(d) 





(3.55) 


Substituting Equation (3.45) in (3.55) yields 


PB@ 





Dhrin =D sy =5.A7 





athe } O(JaR-1,)- G.56) 


Compare Equation (3.56) with (3.47) and notice that the argument of the Q- function is 
slightly smaller now, which results in a tighter upper bound. Substituting Equation (3.56) 


in (3.54) and recalling that R, = 1/2, dice =4 and k = 4 for the dua-4 code, obtains 


P. - Lowel Gar. } (3.57) 


Using the formula in the Appendix for B(d) of dual4 code leads to B(d) = 15, 60, 570, 


2820, and 18165, for d= 4, 5, 6, 7 and 8, respectively. Approximating Equation (3.57) 
with its five first terms, the final formula for P, with P.o, = 95.5% in Figure 3.5 is 


obtained as 


n= [iso i +000 J25%)+5700( 7) 
+28200(/35y,)+181650(./4y, }]. 


Notice that the simple assumption of D‘,, =5A” made the upper bound for P, 


(3.58) 


tighter in 1 dB, which can make a difference when designing a communication link. As 
stated previously, this assumption is valid for 99.5% of the transmissions. However, the 
remaining 0.5% will still follow the original upper bound, which could be considered bad 
if 100% of the data were submitted to it, but considering that only 0.5% of the data will 


ever truly be submitted, this original upper bound can become very acceptable. 
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Figure 3.5. Difference on Performances when Changing D-,, on the Calculation of 


the Upper Bound of P,(d) from D*,, =4A* to D‘,, =5A° for a 16-QAM Channel with 
Dual4 Convolutional Code and Soft-Decision Viterbi Decoding. 


How far can the process of increasing D*., go? Can it be increased to a point 
where 95% of the channel data will follow the new correspondent upper bound? 
Probably yes, but the precise answer to these questions will only depend on the 
application intended for the communication link. Some practical experiments should be 
done in order to define the acceptance threshold. Anyway, this approach provides a 
tremendous enhancement to the upper bound, allowing a more realistic view of the 
behavior of some good, simple and cheap configurations that probably would be 


discarded at first glance if considering just the worst-case upper bound. 
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Recall that Equation (3.52) was created assuming that the second smallest 
possible value for D%,, was chosen. However, for general choices that could be greater 


than the second smallest, this formula will change to 


1 d 
P..(Diy,) = 1 Sire] Py —C, ll? =Diy; ; (3.59) 
i=] l=1 


where P. (Diy) represents the percentage of the total d-symbol combinations that the n- 


th smallest possible value of D%,, will cover. Note that Digs = D?.. by the previous 


symbology. 
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Figure 3.6. | Dpdfof Dj, for a 64-QAM Rectangular Constellation, d = 3, and 
Conservative Choice of Co. 


Figures 3.6 and 3.7 describe two dpdf’s for a same 64-QAM constellation in 


which the first plot corresponds to a path presenting d = 3, while, in the latter, d = 4. The 
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previously mentioned effect of diluting the dpdf of D%,,/ A’ as d increases is easily 
identifiable in these figures. The number of possible values for D‘,, found in both cases 
were 79 and 110, respectively. Notice that in Figure 3.7, since the D{,, domain is greater, 
the overall height of the dpdf is smaller than that in Figure 3.6. The probability of finding 
one value for D%,, that is smaller than a certain fixed threshold reduces when the dpdf 


shortens. 


As previously exposed, since the D%,, determined for a path with d = djiee will be 
even better for another path with d > dyee, Figures 3.6 and 3.7 can be used to study the 
association of a 64-QAM modulator with different convolutional codes presenting 
djree = 3 (Figure 3.6) and djree = 4 (Figure 3.7). For instance, if D*,,/ A* =20 is chosen for 
improving P,(d), by entering the data of the dpdf’s in Equation (3.59), the probabilities 


that the improved upper bound covers all data transmission would be 95.0% for dyree = 3, 
and 97.2% for dyree = 4. Notice that for the first case, in which dyree = 3, there are 12 
values in the horizontal axis of the dpdf in Figure 3.6 that precedes 20, whereas for 


djree = 4, this number increases to 16. 
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Figure 3.7. Dpdfof Dj, for a 64-QAM Rectangular Constellation, d= 4, and 
Conservative Choice of Co. 


Although there are 16 values for D%,, / A’ smaller than 20 to be considered in 
Figure 3.7, these values are very small due to the large possibilities of combining 63 
symbols in 4 slots (63 = 15,752,961). Figure 3.8 shows a zoonrin on the first values of 


the dpdf in Figure 3.7. Notice how small the 16 values are. The sum of all of them is 
0.028, and when applied to Equation (3.59), can result in 97.2%. If, instead of 97.2%, the 


specification for the upper bound was 99.5% of the data, a value of D%,, =13A° would be 


appropriate when looking up the dpdf in Figure 3.8. 
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Figure 3.8. Initial Part of the dpdf of D%,, for a 64-QAM Rectangular Constellation, 
d= 4, and Conservative Choice of Co. 


As a final comment on this topic, this method of improving the upper bound 
for P,(d) results in breaking the complete independence between B (d) and P,(d). 
Although most of the independence remains, when considerations are made about dyce 
applied to P,(d), they create an implicit link between the convolutional code and the 
demodulation scheme. In conclusion, the methodology herein exposed is directly 
dependent on the specific communication system being considered and not solely on an 


isolated P,(d). 
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D. PERFORMANCE ANALYSIS 

Finding good convolutional codes is not easy. They are normally found by 
computational researches on a try-and-error basis. The goal of this Section is not to find 
the best convolutional code that can be applied to a certain M-QAM modulator, but to 
demonstrate how the concepts developed thus far in this thesis can be applied to measure 
realistic performances in communication systems when combining a certain M-QAM 
modulator with a certain convolutional code. It is much more a question of showing 
analysis tools and overall behaviors than obtaining the best results. A straightforward 
way to do that is to run performance analyses over a limited number of interesting target 


systems. 


The focus will be concentrated on three convolutional codes and four different M- 
QAM modulators that will be presented in the following order for didactic purposes: 16- 
QAM, 64-QAM, BPSK, and QPSK. The codes were chosen due to either being well 
known or by having demonstrated outstanding performance during computational 
research. The plots of P, versus E;/No were obtained from simulations taking into account 
the methods for determining B (d) and the formulas for finding the upper bound of P, (d) 
and P) presented thus far in this chapter. The improvements in the upper bound of P, (d), 
the subject of the latter Section, were also part of the results. Although the BPSK channel 
is not a norbinary channel, it has been included here for the purpose of comparison since 
its characteristics are usually well known to the reader and very well cited in the 


literature. 


Sometimes the uncoded 8-QAM performance will be necessary for comparison 
between the options of sending coded or uncoded symbols through the channel. The 8- 
QAM modulator will then be assumed to have the constellation in Figure 3.9, which leads 


to the following expression for its uncoded performance: 


=20| [Prin | _ 2 
no 52(V¥.)- (3.60) 
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Figure 3.9. The 8-QAM Rectangular Constellation 


The three convolutional codes are the following: 


(2k, k) dual-k that is appropriate to work with BPSK, QPSK, 16-QAM and 64-QAM, 
or in general, with all 27*- QAM, and has its B (d) determined by the formula 
developed in the Appendix 


(2, 1) IEEE 802.11a code that can be applied to BPSK, QPSK, 16-QAM and 64- 
QAM, widely employed in wireless networks. The correspondent B (d)’s for each 
modulation scheme was determined by the numerical method discussed in Section 
ILD.2, and shown in Table 2.1 


(4, 3) CHL code from Ref. [7]. This code will be referenced from this point forward 
as Chang-Hwang-Lin, or CHL for short, in honor of the three discoverers. It is 
appropriate for 16-QAM only and counts on having, as an advantage, a relatively 
high code rate compared with the other two (2, 1) codes. By this latter characteristic, 
CHL does not compromise the maximum bit rate capacity available in the channel 
very much. Table 3.1 depicts the main characteristics of this code, where G is the 
octal format representation of the transfer- function matrix G(D) as defined in Chapter 
IL. 
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Table 3.1. Main Characteristics of the (4, 3) Chang-Hwang-Lin Code Employed in 
this Thesis. 


1. 16-QAM with Dual-4, IEEE 802.11a and CHL 


Figure 3.10 shows the various performances for a 16-QAM with a dual-4 encoder. 
The performance of an uncoded QPSK channel was also plotted because it represents an 
alternative to sending the same 2 information bits contained in the 4-bit 16-QAM coded 
symbol during one transmission. Recall that R, = 1/2 for duakk codes. Keeping this 
alternative in mind, the first goal of a good convolutional code with R, = 1/2 applied to a 


16-QAM channel is to exceed the performance of the uncoded QPSK. 


Looking at Figure 3.10, the dual4 is able to exceed the uncoded QPSK during 
94.4% to 97.5% of all transmissions for values of 2, greater than 8 dB. Notice that the 


consideration of only the first upper bound corresponding to Peo, = 100% would at once 


discard the dual4 code. 
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Figure 3.10. 16-QAM with (8, 4) Dual 4 and Viterbi Soft-Decision Decoding and 
Improvements in the Upper Bound. 


Another interesting behavior to note is how fast the upper bound tightens itself at 
the first improvement steps. Although the amount of reduction on P.o, increases during 
the last improvements, the corresponding variation in the upper bound becomes smaller. 
This behavior can be explained by looking at the dpdf of the 16-QAM constellation for 
djree = 4 in Figure 3.4. The first point corresponds to a small probability that makes P,o, 
vary just a few small steps. On the other hand, the subsequent improvements reflect at 
least in 3% of the changes in P.o,. The analysis of the dpdf determines what can be 


expected from the improvements. 


Table 3.2 shows the correspondence between the values of D?,, and Peoy 
considered for each improve ment in P,. These values can be achieved using the dpdf in 
Figure 3.4 and the method described in Section IIL-C.3. The first four were used in the 


plot in Figure 3.10. 


aD 





fare (%) 


Table 3.2. Values of D%,, / A’ Correspondent to Different Probabilities of Coverage 
(Pcov) for a 16-QAM Rectangular Constellation Using Convolutional Code with dye = 4. 


Finally, notice that due to improvements, it is possible to tighten the original 


upper bound up to 2.5 dB and still have a reliable curve that covers 94.4 % of all data in 
the channel. 

Figure 3.11 illustrates the behavior of the (2, 1) IEEE 802.1 1a code when applied 
to a 16-QAM channel. In general, the behavior is very close to the duak4, being slightly 


worse. 
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Figure 3.11. 16-QAM with (2, 1) IEEE 802.11a Convolutional Code, Viterbi Soft- 
Decision Decoding, and Improvements in the Upper Bound. 
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All the considerations made for the dual 4 are still valid except that only one bit is shifted 


per time into the (2, 1) IEEE 802.11a encoder contrasting with the four bits of the dual-4. 


Figure 3.12 depicts the behavior of the (4, 3) CHL code working with the same 
16-QAM channel. The performance now is almost 1.5 dB better than the other two cases. 
Considering that this code consumes only 25% of the channel bandwidth in contrast to 
the 50% achieved by the other two, the obtained result becomes much better. In fact, 
notice that the basis for the performance comparison in Figure 3.12 is the uncoded 8- 
QAM curve instead of the uncoded 4-QAM used on both (8, 4) dual4 and (2, 1) IEEE 
802.11a. The difference occurs because now there are three information bits per 
transmitted coded symbol instead of two, which could be done through an uncoded 8- 


QAM channel. 


The relative step improvements in the upper bound when applying the (4, 3) CHL 
can be seen to be the same in Figure 3.12 as those obtained in Figures 3.10 and 3.11. As 
was demonstrated previously, the effect of the convolutional code on those improvements 
depends only on its dye. Since, coincidently, all three codes exposed here present dyree = 4 
when working with the 16-QAM channel, the improvement steps will be exactly the same 


for all of them. 


Clearly, the best choice among all three codes to associate it with a 16-QAM 
modulator falls to the (4, 3) CHL. When trying to express this performance in one single 
metric, note in Figure 3.12 that only 8 dB are required for the SNR to achieve P, = 10° 
for 94.4% of the data. It represents a gain of 4.5 dB when compared with the uncoded 8- 
QAM at the expense of only 25% of bandwidth loss. 


Another remarkable fact about (4, 3) CHL code is its K = 4, considered a relative 
small constraint length for good codes. Just for the sake of comparison, IEEE 802.11a 
code presents K = 7, although its shift registers size is equal to 1, which is smaller than 3 


used by the (4, 3) CHL. 
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Figure 3.12. 16-QAM with (4, 3) Chang-Hwang-Lin Convolutional Code, Viterbi Soft- 
Decision Decoding, and Improvements in the Upper Bound. 


2. 64-QAM with Dual-6 and IEEE 802.11a 


Figure 3.13 illustrates some performances of a 64-QAM channel associated with a 
dual6 encoder. In this case, since R, = 1/2 and g = 6 bits, the alternative uncoded channel 


used for comparison is the 8-QAM. 


Just looking at the performance of the upper bound that covers 100% of the data 
in the channel clearly indicates that the dua-6 code does not lead to significant 
enhancements over the uncoded 64-QAM channel. On the contrary, the small gain it 
provides is much less than the choice of an uncoded 8-QAM. However, considering the 
improved upper bounds, this apparently weak dual-6 code becomes a very interesting 
option. Notice especially the first improved upper bound that covers 99.6% of the data. It 
represents about 5 dB in the coding gain analysis along a range that practically goes 
through all possible sequences in the channel. In fact, considering only the curve that 


Poy = 99.6%, it overcomes most of the gain obtained if choosing the uncoded 8-QAM. 
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The characteristic of the large shift in the upper bound for 64-QAM when 
compared with 16-QAM is very comprehensible when analyzing the dpdf of D4, for 


both constellations shown in Figures 3.4 and 3.7, respectively. As explained in Section 
Il.C.3, when commenting about Figure 3.7, the greater granularity presented by the dpdf 


of 64-QAM points out toa D%,, for covering 99.6 % of all possible sequences much 
farther from 4A’ (D%,, for 100%) than this same parameter for 16-QAM. In other words, 
the D%,, that makes P.., = 99.6% is much greater for the 64-QAM than for the 16-QAM 
constellation, even though both present 4A” as their minimum possible value for D’,, . 
The farther the considered D%,, is from the minimum, the larger the first improved upper 


bound will shift from the original one. 
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Figure 3.13. 64-QAM with (12, 6) Dua-+6 Convolutional Code, Viterbi Soft-Decision 
Decoding, and Improvements in the Upper Bound. 
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min ? 


Table 3.3 shows the correspondence between values of D%,, chosen to be D 


and Po, for each improvement on Py). 
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Table 3.3. Values of D%,, /A* Correspondent to Different Probabilities of Coverage 
(Poy) for a 64-QAM Rectangular Constellation Using Convolutional Code with dye = 4. 


Finally, it is remarkable how the dual 6 code demonstrates an excellent 
performance with 64-QAM after all improvements. Notice on the curve for Peo, = 95% 
that it requires only 8.5 dB to provide a bit error rate of 10°. This performance is 
comparable with that shown for 16-QAM with (4, 3) CHL code in the latter Section. 
Furthermore, this result is enhanced by the fact that the dual6 could be categorized as a 
poor code to work with 64-QAM when analyzing only the most conservative upper 


bound, which would lead to a complete distortion of a realistic behavior. 


Figure 3.14 shows the behavior of the (2, 1) IEEE 802.11a code when applied to a 
64-QAM channel. The overall performance of this code can be seen to be approximately 


2 dB worse than the dual6 when considering the improved upper bounds. 


The same comments related to the shifts on the improved upper bounds made for 
duaL6 is also applicable for the (2, 1) IEEE 802.11a. However, notice that djee now is 3 
instead of 4. Consequently, the dpdf for this case is the one shown in Figure 3.6, which 
presents less granularity than that in Figure 3.7 for djee = 4. Less granularity means 
smaller distances from 4A? for the D7‘ telated with the first improved upper bound, 
which leads to a smaller shift. In fact, notice that the first shift in the upper bound for the 
(2, 1) IEEE 802.11a is 4 dB and reaches P.o, = 99.4%, while approximately this same 


parameter for the dual6 code reaches 5 dB. 


The small djee presented by the (2, 1) IEEE 802.11a when associated with 64- 
QAM can be pointed out to be the cause of this code being outperformed by the duak6. 
However, notice that for Peo, = 95% the (2, 1) IEEE 802.11a code does not disappoint, 
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achieving 8 dB of coding gain from the uncoded 16-QAM curve at P, = 10°, and 4.5 dB 
when compared with the uncoded 8-QAM. 
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Figure 3.14. 64-QAM with (2, 1) IEEE 802.11a Convolutional Code, Viterbi Soft- 
Decision Decoding, and Improvements in the Upper Bound. 


Figure 3.15 again demonstrates the performance of the two previously discussed 
codes but are no w plotted together and only for the most conservative upper bound. 
Clearly, it can be seen that the dual 6 is better than the IEEE 802.11a. Moreover, a 
performance analysis based only on these plots would lead to the discarding of both 
codes in favor of the uncoded 8-QAM option, which was seen to be a wrong conclusion. 
This last fact highlights the great usefulness of the upper bound in improving 
methodology for large constellations. The larger the constellation and the dy-ce, the better 


will be the improved upper bounds. 
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Figure 3.15. 64-QAM with (12, 6) Duak6 and (2, 1) IEEE 802.11a Convolutional 
Codes for Peo, = 100% Using Viterbi Soft-Decision Decoding. 


3. BPSK with Dual-1 and IEEE 802.11a 


Figure 3.16 reports the performance of the (2, 1) duak1 and (2, 1) IEEE 802.11a 
codes in a BPSK channel. In this case, improvements in the upper bound are no longer 
possible due to the nonexistence of a possible dpdf of D‘,, for BPSK. Recall that BPSK 
presents only two symbols and, regardless of the symbol chosen to be Co, there is only 
= 4A’, 


one option left for a possible erroneous coded symbol, which leads to D‘,, = D- 


for every dyree. 

Clearly, it can be seen that the IEEE 802.11a is much better than the duak6 in 
providing almost twice the coding gain in dB. Notice that the relative performance here is 
the opposite of the 16-QAM and 64-QAM cases for the same analysis. Recall from Table 
2.1 that dyree for the IEEE 802.11a decreases as M increases, obtaining the maximum 


value of 10 exactly for BPSK. Whereas, for the duakk code, since the number of memory 
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elements in the shift-registers accompany the increasing of M, its dye remains constant at 
4, When dyree is higher for the IEEE 802.11a code than for the dua-k, the first code 


outperforms the second. 
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Figure 3.16. BPSK with (2, 1) Duak1 and (2, 1) IEEE 802.11a Convolutional Codes 
Using Viterbi Soft-Decision Decoding. 


4. QPSK with Dual-2 and IEEE 802.11a 


The QPSK modulation is midway between BPSK and 16-QAM in terms of the 
possibility of improvements. Contrary to BPSK, it can present a dpdf of D%,, but is not as 
prolific as a 16 or 64-QAM. The QPSK constellation is very symmetrical, which allows 
any choice of the Co symbol to have the same result. In other words, there is no 
conservative choice for Co in the QPSK rectangular constellation. Also, once Co is 
chosen, among the other three remaining symbols, only one of them presents an 
Euclidean distance different than the minimum, which allows the existence of a timid 


dpdf that can only obtain higher granularities if d becomes too large. Figure 3.17 depicts 
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the dpdf’s for the values of djee equals to 4 and 6, corresponding to the dual 2 and IEEE 
802.11a codes, respectively. 


The dpdf for d = 4 can be seen to have very small granularity presenting only five 


2 
AV? 


points. Therefore, the minimum value for D>, , 4A’, can be present in 20% of the possible 


sequences in the channel. This means that the minimum improvement achievable for the 
upper bound will embrace only 80% of the possible transmission, which can considerably 
damage theoretical performance forecasts. On the other hand, when d increases just a 
little bit and becomes six, the dpdf granularity increases to seven points and the minimum 
improvement can now cover 91.2% of the data. Notice in Figure 3.17 that there is 8.8% 
of probability of happening the minimum average distance for d = 6. Although, 

Peoy = 91.2% is much less than the minimum of 95% considered before for 16-QAM and 
64-QAM,, it is possible to imagine a design considering this value if it yields considerable 
improvement in the upper bound. The restrictions discussed in this paragraph clearly 
demonstrate how the efficiency of the upper bound improvement directly depends on the 


magnitude of M. 
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Figure 3.17. Dpdf of D%,, for a QPSK Rectangular Constellation Considering Two 
Values of d, 4 and 6. 
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As a result of the previous issues about QPSK dpdf’s, the performance analysis 
does not consider improvements for the upper bound of the dual2 code because it 
presents djrce = 4. The only improvement to consider is for the IEEE 802.11a code, which 
shows dyree = 6 when working with QPSK, as can be seen in Table 2.1. Actually, this 


procedure is limited to a first step in which D‘,, = 4.6667 A* that corresponds to 


Peoy = 91.2%. Figure 3.18 illustrates the QPSK performance analysis for the two 
considered codes. Since R, = 1/2 for both codes, the alternative uncoded channel that 
QPSK can compare with is the BPSK, which can be seen to perform the same as the 
uncoded QPSK. In fact, QPSK works like BPSK operating in two completely 
independent orthogonal channels, which explains the identity of these metrics. This 
means, in practice, that there is no alternative uncoded channel to use instead of coded 
QPSK. Consequently, any possible coding gain will be very welcome. The reader can see 
how to this point the QPSK case differs from the 16-QAM and 64-QAM, presenting at 


the same time, characteristics of a binary and a nombinary modulation technique. 


Looking at Figure 3.18 for P.o, = 100%, the IEEE 802.11a code can be seen to 
provide more coding gain than the dual-2 code, even though the difference between these 
gains is now smaller than for the BPSK case. Notice that the djc¢e for the IEEE 802.1 1a 
code drops from ten to six when the modulation changes from BPSK to QPSK. However, 
six is still greater than four, djee for the duak2 code, pushing IEEE 802.1 1a to overcome 
the performance of dual2 for QPSK. 


Concerning the improvement in the upper bound for the IEEE 802.11a code in 
Figure 3.18, the resultant shift can be seen to be very small, almost 0.7 dB. This is 
justified because D%,, (= 4.6667A’, for this case) is too close to the minimum distance 
4A’. The consequent loss in P.o, does not compensate for this small improvement and the 
best decision should be to disregard this curve. On the other hand, for the performance of 
TEEE 802.11a with P,., = 100%, the coding gain at P, = 10°° can be seen to be 4 dB, 
which is less than the 5.5 dB obtained for BPSK, but still a good coding gain. Notice that 
the dual-2 code provides only 2.5 dB for this gain. 
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Figure 3.18. QPSK with (4, 2) Duak2 and (2, 1) IEEE 802.11a Convolutional Codes 
Using Viterbi Soft-Decision Decoding. 


E. SUMMARY 

Chapter II dealt only with the pure AWGN scenario, deriving the upper bound of 
P,(d) for any M-ary coherent demodulator followed by a SDD. This upper bound was 
customized for M-QAM and generated a model for determining P, for any M-QAM 
modulation with SDD. The chapter also demonstrated how to take advantage of the large 
number of symbols existing in M-QAM and improve P, in an attempt to tighten it to 
more realistic data in the channel. As a practic al application, the performance analyses of 
16-QAM, 64-QAM, BPSK and QPSK, associated with up to three different convolutional 
codes were shown. Chapter IV will illustrate how to adapt all results obtained in this 


chapter to derive P, for Nakagami fading channels. 
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IV. PERFORMANCE ANALYSIS IN A NAKAGAMI-M CHANNEL 


All issues discussed in Chapter III only account for one type of random variable 
in the receiving process, which is related to the AWGN. Although this approach is 
suitable for many types of communication links, such as free space and wired 
communications, unfortunately it does not completely satisfy the requirements of some of 
the most common applications. Whenever obstacles are present between the transmitter 
and the receiver that can create scattering and multipath, the communication engineer 
must consider the effect of more random variables in the channel model. This is the case 
with wireless networks and wireless telephony, in which walls, buildings, cars, trees and 
all sort of obstacles can be between the two communicating points. Other examples are 
the radio ionospheric communications in the HF band and ionospheric forward scatter in 
VHF. In these cases, the ionosphere is responsible for reflecting and scattering the 
signal, and allowing it to reach the receiver when there is no line-of-sight (LOS). This 
chapter makes a parallel distinction between the pure AWGN scenario studied in Chapter 
III and the scenario in which all the effects of midway elements can be modeled by a 


Nakagami-m fading channel. 


A. NAKAGAMI-M CHANNEL 

Many elements can affect signal propagation, such as multipath, Doppler effects, 
shadowing by surrounding objects and changes in environmental characteristics. The 
final result is a channel that has a randomly time variant impulse response. In other 
words, if the same impulse is transmitted successive times, the correspondent received 
signals will be different each time. The differences can be in the number, in the relative 
delays and in the shape of the received signals related to each transmitted impulse. Such 
channels are named time-variant channels and can be categorized in terms of the 


frequency selectivity and fading characteristics they introduce in the signal. 


The best model for a certain communication link depends on its bandwidth and 
the environment in which the communication end-points will operate. A large number of 


channels present the characteristic of non-distorting signal harmonics, i.e., the frequency 
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response of the channel is flat for the entire bandwidth. Also, the fading conditions in 
those channels vary so slowly during a symbol interval (7s) that the resultant attenuation 
can be considered constant along each Ts. Even being constant during a certain Ts, notice 
that the fading conditions can vary from symbol to symbol. Those channels form a 
category named frequency-nonselective slowly fading channels (or flat-slow-fading 
channels) and result in a multiplicative distortion of the transmitted signal by a random 


variable attenuation factor O& as in Ref. [3]. 


The Nakagami-m channel to be considered in this thesis is a frequency nor 


selective slowly fading channel in which o is characterized statistically by a Nakagami- 


m distribution and has E [a | =1. It is the best statistical model for this category of 


channels and is applicable in many communication systems that work under fading 
effects, such as wireless networks and wireless telephony. A Nakagami-m channel is the 
best fit for data signals received in urban radio multipath channels [3], but it also 
embraces other models in special cases such as Rayleigh, utilized when there is no LOS, 
and Ricean distributed, appropriate in the presence of LOS. Notice that the channel noise 
remains the same AWGN as previously seen. Therefore, in the subsequent analyses, two 
types of random variables will be considered,  , that follows the Nakagami-m 
distribution, and n,;, for i =1 and 2, that follows the Gaussian distribution and represents 
the two components of the AWGN that are present in the two orthogonal channels of the 


coherent M-QAM demodulator. 


Applying the Nakagami-m distribution [3] to leads to 


im —mo" 


Z m (2m-1) [a” | 


“Fo Eo 





F,,(@) (4.1) 


where I'(m) is the Gamma function, which is defined as 


I(m) = er ler ae m>0O. 
0 
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The parameter m 1s called fading figure and can be seen as the ratio 


m= ee ie (4.2) 


E| (a?-E [a })| 2 
Notice that the denominator of m is the variance of the second moment of @. 


Recall that 
Ela? ]=1. (4.3) 
Substituting Equation (4.3) in Equations (4.1) and (4.2) yields 


(m)" Derm 


fy (@) (4.4) 


~ Fan) 
and 


pees. ype (4.5) 


B[ («*-1)| 2. 


Notice in Equation (4.4) that when m=1, f, (a) reduces to 20 e° , which is a 


Rayleigh distribution. Also, when a” remains constant and equal to 1 for every received 
symbol, no variance exists in the attenuation factor, which causes Equation (4.5) to obtain 


m =o, This is what characterizes a no fading scenario in which only AWGN affects the 


receiving process. All the plots obtained in Chapter III are included in this case. 


As stated previously, the attenuation factor & imposes a multiplicative distortion 
in the transmitted signal in the cases of frequency- nonselective slowly fading channels. 


Thus, the new received signal can be represented as 
r(t) =Q s(t) + n(t), (4.6) 
where s(t) is the transmitted signal component, and n(t) characterizes the AWGN. 


By Equation (4.6), the average signal energy contained in one received symbol in 
the presence of attenuation factor & can be expressed as 
S185 (4.7) 
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where €, is the average signal energy contained in one received symbol if no attenuation 


exists in the channel. 


Since the noise remains the same, Equation (4.7) leads to 
Yon =O°Vs, (4.8) 


where y,,, is the average SNR in a symbol for a given and y, is the average SNR ina 


symbol when there is no attenuation. 


Also, the overall average SNR in a symbol can be expressed as 
To = Ela? |y,=Y,5- (4.9) 
Substituting Equation (4.9) in (4.8), results in 
You =O Fy « (4.10) 


Since o& 2 0, Equation (4.10) leads to a one-to-one mapping between and y,, . 


Hence, the pdf of y,,, can be straightforwardly obtained by 


Src Lin} (2), ey: (4.11) 


lay, /da] 


Working with Equation (4.10), the following relationships stand 


o = [oe (4.12) 
VT 


oe 


and 


= 207... 4.13 
ef van (4.13) 





Applying Equations (4.4), (4.12) and (4.13) into (4.11) leads to 


m TINY S01 


m 
i. 03) ere © (4.14) 
ee Hoya) 
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The interpretation of the latter equation can be as follows. If a sequence of a 


same symbol is sent in a Nakagami-m channel with fading figure m, the set of all SNR of 
each of the equivalent received symbols, y,,,, would constitute a random variable that 
varies accordingly to the pdf provided by Equation (4.14). Notice that the reason in the 
previous description for a same symbol being sent is only to simplify the understanding. 
Recall that in M-QAM modulation, different symbols can have different amplitudes 
leading to different SNR. This explains the application of the term “average” for the 
definition of Y,;,, in Equation (4.8). Also note that even though y,,, is a random variable, 
Ysq n Equation (4.14) is a deterministic value equivalent to the SNR that those symbols 


would obtain in the absence of attenuation in the channel. 


B. PROBABILITY OF BIT ERROR FOR UNCODED M-QAM IN 
NAKAGAMI-M CHANNEL 


The performance analyses presented in this thesis always compare the probability 
of bit error P, of aM-QAM modulation that makes use of a certain convolutional FEC 
versus the option of sending uncoded symbols through the same channel. Thus, in order 
to study the M-QAM performances in a Nakagami-m channel, it is necessary to obtain 


the corresponding expressions relating P, with the average SNR per information bit, 


called y,,. Recall that in the case of Nakagami-m channels, since each symbol can be 
received with different attenuations, the SNR per information bit, y,, is a random 
variable that can be expressed as a function of y,,,. Therefore, in order to compare 
performances in Nakagami-m channels with those previously obtained, the average Y, 
must be used since it is a deterministic value and represents the same y, employed for 


the pure AWGN case. The difference in the notations will be kept to emphasize the 


averaging reference. 


In a pure AWGN scenario, the following expressions stand for determining the 


upper bound of P, of an optimum receiver for a M-QAM uncoded signal [3]: 


a 


1 
PS =P.. (4.15) 


and 





Py = vol Prin } (4.16) 


where gq is the number of bits per symbol, Py is the probability of symbol error, N,, is the 


largest number of neighboring points that are at distance D-,, 


from any constellation 


2 


points, and D-,, is the minimum Euclidean distance between signal points in the 


constellation. 
Equation (4.16) can be rewritten as 


De é D2, 
P,=N Ain “8 |= N miny |, 4.17 
M Z| 2 : | Z| “| ( ) 





where €, is the average energy of one received coded symbol. 


Consider a Nakagami-m channel that presented a certain attenuation factor & 
applied to a received symbol and correctly compensated for by the detector, i.e., a 
specific given 1. Due to the energy issues discussed in Equation (4.7), Py can now be 
computed for this specific symbol as a function of & by means of 

D2 


P,(a)=Nuy | Se Us } (4.18) 


By substituting Equation (4.8) in (4.18), Py can be expressed as a function of Y,,, 


as 


D?. 
Py (Ysa) = Nu Q| IY 50. |- (4.19) 
2€, 


2 


Notice that Equation (4.19) is valid only if y,,, is known, i.e., this is a conditional 
Py. In order to determine the average upper bound for Py, an integration with the pdf of 


the random variable Y,, is necessary and is given as 


Py = (2 Cs) Jy (Ve ) AY se . (4.20) 


Substituting Equations (4.19) and (4.14) in (4.20) leads to 








+o0 D?. m” asa 
Pur y Nu | Q = Vso Tm vo eu AY 5 i (4.21) 
2€5 Tn) (Tsu) 

The analytical solution for this integration is not easily obtained. However, 
following the solution given in [8] for a similar integral when analyzing the special case 
of BPSK, after a development involving Gauss’ hypergeometric function and 


Pochhammer’s Symbol, the following result can be achieved, 





———— (4.22) 
1+ Dyin ¥ so 
4esm 


From the previous discussion about Y,, at the beginning of this Section, notice that 


Vs0 = 47> (4.23) 


Substituting Equations (4.22) and (4.23) in (4.15), the final equation for P» is 


obtained as 
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ere 
| 5 | 


) MEP msn + Death , } 


4e «m 


Di Ns 
4e «m 


(4.24) 





Notice that Equation (4.24) is generic and can be applied to every M-QAM 


constellation. 


In Section III.C.2 it was demonstrated that each M-QAM constellation has a 


different ratio D°,,/€, that depends on M. Also, for rectangular constellations in which 
M =2‘ and q is even, there is a special method in [3] for precisely determining Ny that 


leads to 


Ni (4.25) 


For odd qg, Ny approaches the average number of neighboring points that are at a 


distance D? 


min * 


Table 4.1 summarizes the ratios D*,, /€, as well as the corresponding 


Nw’s for each M-QAM rectangular constellation considered in this thesis. The only 
modulator that presents odd q in Table 4.1 is the 8-QAM, which its constellation is 


depicted in Figure 3.9. 
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Table 4.1. Useful Coefficients for Rectangular M-QAM Constellations. 


For the special cases of rectangular constellations, the values in Table 4.1 should 
be used into Equation (4.24) to obtain the various performances. For BPSK and QPSK, 
since N,,/q =D*. q/(4e€,) =1 for both, Equation (4.24) provides exactly the same 


expressions for these two modulations, which matches the equivalent behavior in the 


AWGN only scenario. 


Figures 4.1, 4.2 and 4.3 show the P» for the uncoded BPSK/QPSK, 16-QAM and 
64-QAM, respectively. When m= 1, as discussed previously, the Rayleigh distribution 
for the attenuation is observed. Recall that small values of m means high variances of 
a”, which signifies more severe conditions of fading. Note the presence of a tremendous 
drop in the performance for small m’s (0.5, 1 and 1.5) in all three plots. Moreover, in this 
high fading situation, there is almost no difference in performances among BPSK/QPSK, 


16-QAM and 64-QAM. They become more evident for larger m’s. 
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Figure 4.1. | Uncoded BPSK/QPSK in Nakagami-m Channel. 
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Figure 4.2. | Uncoded 16-QAM in Nakagami-m Channel. 
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Figure 4.3. | Uncoded 64-QAM in Nakagami-m Channel. 


C. PROBABILITY OF BIT ERROR FOR M-QAM WITH SDD IN 
NAKAGAMI-M CHANNEL 


The principles of determining the upper bound of the probability of bit error P, 
for M-QAM when coding is being used in a Nakagami-m channel are exactly the same as 
previously seen for AWGN, indicating that the general equation shown below remains 


valid, 
i 
P,Sz y" B(d) Pd). (4.26) 
d= free 


The parameter B (d), which corresponds to an exclusive characteristic of the code 
being utilized, is the same seen in Chapter II since the code itself did not change. The 
only parameter that varies is P, (d). As seen in the latter Section, P, of an uncoded 
symbol increases tremendously for severe fading conditions (m < 1.5). Therefore, the 


occurrence of a larger number of wrong paths should be naturally expected through the 


decoding process. 
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In this development, the same thread of thoughts followed in Section II.C.1 will 
be adapted to the Nakagami conditions in order to obtain the expressions for P, (d). The 
final P;, will be obtained by substituting this new P, (d) and B (d) , already computed in 


Chapter II, in Equation (4.26). 


1. Upper Bound of P2(d) for any M-ary Coherent Demodulator and SDD 
in Nakagami-m Channel 


Although the channel is now a frequency nonselective and slowly fading 
Nakagami channel, the receiver is still the same and is shown in Figure 3.2. Consider the 
reception of a sequence of coded symbols that presented various attenuation factors o's 
correctly compensated for by the detector, 1.e., consider a given sequence of a's . Each 
Q is related to only one received symbol in the sequence and is considered constant 


during the entire symbol period Ts. Therefore, the notation © ,,, will refer to a particular 
©. that affected the reception of the m™ symbol of the j"" branch in the trellis. 


The expression for the received signal r(t) changes to the one expressed in 


Equation (4.6), which can be written in a vector form as 


Ry =O pCi +N py (4.27) 


jm jm ~ jm 


The branch metric yw” , associated with the j‘” branch of the i” path in the trellis, 


also changes slightly. Since the decoder should compensate for the various © ,,, in the 
signal during the decoding process, this new metric becomes 
log 7 .Re Rs Racks Ry lOue, OC. 026% 080 ))) 
Hi; 0g R jl ONjp2 SS73 °° Th IN ji~ Jl j2— j2 p37 73 ost JN ~ jN 
log fC Re mal ZN NOC mal 2 aN} 


WP Sloe fe Ry Mi oyaCy. 1): (4.28) 


Following the same development as in Section II.C.1, Equation (4.28) becomes 
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N 
w= K,+K, VR, —O ,C5p I? (4.29) 


Q in jm 
m=l 


where 








1 
K, =Nloa{ ===> and K,=- 


Considering the correlation metric adopted by the SDD, the probability of 
choosing a wrong path in the trellis in a pairwise comparison of the all- zero path (correct 


path) with a wrong path 1, which presents distance d from the all-zero path, given a 


certain set of o,,,, is represented by P, (d ) and can now be written as 


P,(d,{o,,}) =Pr[CM” >CM® | 
P,(d,{o.,,,}) =Pr yy (lz, =O, CO IP=WR, — OG COW?) <0]. 4.30) 
j=l m= 


Considering again that a sequence of all- zero symbols is the one being 


transmitted, hence 
Ci, =Cy, Vj andm, (4.31) 
and Equation (4.27) becomes 


R in =% imCo +N im (4.32) 


jm 


Substituting Equations (4.31) and (4.32) in (4.30) yields 


A (af) =P EY (bao 4N, SCO IPS Wey. Ns 2G, I")so| 


j=l m=i 


j=l m= 


P,(d,{o.,,,}) = rfEX( Ob jp (CO —Cy) —N jy WIN, HI?) <o| 


ply 


j=l m= 


Pale} =PH EP ? Wc® -C, 1? -20,,(€% - G).Nn) $0] (4.33) 
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Recalling that C” == =C, , for all combinations of j and m except for d 


jm 


symbols, Equation (4.33) can be rewritten as follows 


l=1 


d 
P,(d,{a.,}) = Pe] Solar IC® —C, I? -20,(C®- C,)N,)s (4.34) 


Expanding the noise vector in its orthogonal components in Equation (4.34) 


results in 


d 
P,(d,{o,}) =P Sor IC® —C, W? —204, CP -c,,) 1, -20,(c9 —¢,) n,}<0 (4.35) 


=1 


Notice that the set of all @,, for / =1 to d, is considered “given” by the first 


assumption of this development. Furthermore, the sequence of C ae for /=1 to d, that 
compounds path 1, is also known by the decoder. Therefore, due to nj and np being 
independent- identical distributed zero-mean Gaussian r.v.'s and the other terms being 
deterministic values, the left side of the inequality in Equation (4.35) can be seen to be a 


Gaussian random variable X expressed as 


xX =" (ER, Yo! Oy } (4.36) 
where 


xX, =a 1C°-C,, (4.37) 
and 


62, = 40? (cP 04)? +(cP -€,,)7)o? =40711C® -C, IP 6? 


=4X 0". (4.38) 


(1) 


Here an upper bound can be established by assuming that all C;’’s, for every 1, 


are equal to a worst case symbol, Cw, in such a way that P, (d ) becomes the highest 


possible. 


If Cwis a constant vector for every /, then Equation (4.36) becomes 
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d d 
X=N [i Cy —Cy IP Ya?, 4G, —C, IP o* Pa; } (4.39) 
1=1 [=1 


In order to maximize the probability in Equation (4.35) and obtain an upper bound 
for P, (d) , Cwmust be the closest symbol to Co. The same theorem and arguments from 


Section III.C.1 can be used to demonstrate this fact. This situation leads to 


IIC,, —C, IP = De 


O,min * 


(4.40) 


Substituting Equation (3.26) in (3.24) and the result in Equation (4.35) yields 


1=1 1=1 


d d 
P,(d,{a.,}) = Pr[X <0] = Pi vi. Ler Dino? La }s0] 





D* 


d 
2 
0,min Ya; 





d 
2 2 2 
AD) min 20; 
[=] 


which implies 


, (4.41) 





where the last step stands due to o* =N,/2. 


Notice that Do min can vary, depending on the Co chosen to be the zero-symbol in 
the constellation. Therefore, in order to guarantee that Equation (3.3 1) will be the real 


upper bound for a given set of &.,, Co should be chosen to be the one that yields 
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D, 


0,min 


= Dyin » 


where D,nin iS the minimum Euclidean distance in the constellation. 


Equation (3.31) can be rewritten as 


P, (d,{ar, }) =Q 


(4.42) 





All the definitions made in this chapter related to symbol energy are still valid for 
coded symbol energy. However, the letter “C” will be added in the notations to 


emphasize the coded symbol scenario. The following expressions stands 


Cosa, = Oe es > (4.43) 
Ycsa, = 06;Y cs > (4.44) 
Tes, = E[ 07 |Yes =Ves> (4.45) 
and 
Yeso, = qR-Y, (4.46) 


Equation (3.32) can be rearranged as 


D.. 


P,(d,{a,}) =O Prin cs 5 oy? =Q Pain y Va; 
2€-5 Ny St | eh 


P,(d,{a,}) = of |e mi Se | (4.47) 


Replacing Equation (4.44) in (4.47) yields 





cs 1=1 


P,(d.{Yc,}) = olf? Paw Sin | (4.48) 


Assume a random variable z defined as 


d 
22 Dye 4.49) 
l=1 
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Substituting Equation (4.49) in (4.48) leads to 





nla)=0[ min -} (4.50) 


In order to determine the average upper bound for P, (d ) , an integration with the 
pdf of the random variable z is necessary and is given as 
P(d)=[_ BC, 2) fe(z)ae. (4.51) 
Notice that z is the sum of d independent-identicaldistributed random variables 
Ycsa, > Which have pdf already derived in Section IV.A and equal to 
=I Csoy 


= ¥¢55,°" 2 7 (4.52) 
T(m) (Fes, ) 


m 


Fy su (Yew ) = 


Following the solution in Ref. [8] for a similar development involving BPSK and 
SDD, by working with Laplace transforms, it can be demonstrated that the pdf of Z will 
be 


—mz 
md z 


fiz = m — gil) eto : (4.53) 
Td) (Fesa, ) 





Applying Equation (4.23) in (4.53) yields 


—mz 
md mz 


= i ep che (4.54) 
Td) (GR) 


fz(z) 


Substituting Equations (4.50) and (4.54) in (4.20) leads to 


Ee D?. md = 
P,(d)= [7 O) Jom < \——*>____, cm ee de. (4.55) 
~(N2€ cs JP (md) (GR F,) 
Note how Equation (4.55) resembles Equation (4.21). There is no difference with 
respect to the integration variable. Consequently, the analytical solution can be obtained 


by following the steps in [8], as for Equation (4.21). The final result can be expressed as 
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I) md + — 
P,(d)= + 2 =e 
1+( Pane | MET om 4 Pin 
4E osm 4e—5m 





(4.56) 


[| (ma +1+n) 


n=0 


Due to the orthogonalities existent in the demodulating process, all the 
assumptions thus far in this Section for a bidimensional coherent demodulator are also 
valid for mult+dimensional M-ary modulations such as M-FSK. The unidimensional 
types, such as BPSK, can be viewed as special cases that are also embraced by the 
previous development. Consequently, Equation (4.56) can be considered a general 


formula for the upper bound of P, (d ) that can be applied to every M-ary coherent 


demodulator followed by a soft-decision Viterbi decoder in a Nakagami-m channel. 


2. Upper Bound of P2(d) for M-QAM and SDD in Nakagami-m Channel 


For the special cases of M-QAM rectangular constellations, P, (d ) can be obtained 
by simply substituting the correspondent values of Table 4.1 in Equation (4.56). Notice 


that ““D? 


min 


/e,” in Table 4.1 can be used for D*,,/€ ., in Equation (4.56) since the coded 
symbol is the one that really modulates the carrier, performing the role of one regular 


symbol in uncoded systems. 
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After consulting Table 4.1, the expression D..g /(4€,s ) leads to the following 


values to be directly applied into Equation (4.56) for BPSK, QPSK, 8-QAM, 16-QAM, 
and 64-QAM, respectively: 1, 1, 1/2, 2/5 and 1/7. Notice that BPSK and QPSK provide 


the same result for P, (d ) again. 


3. Improvements on the Upper Bound of P2(d) in Nakagami-m Channel 


All considerations regarding improvements on the upper bound of P, (d) made in 


Section II.C.3 are still valid for a Nakagami-m channel. The previous study was based 
on the nature of the constellation, the probability of transmitting some specific code 


sequences, and djce. There is no relation with the type of channel being used. Therefore, 


2 


instead of using D%,,/€., in Equation (4.56), D‘y /€¢, should be employed. The 


2 


min * 


parameter D‘,, represents the square of the new average distance to be used as D 


A practical procedure for applying a certain improvement represented by 


Di, | A’ can be derived as follows: 


Diy ss PisyPas = Prin Day = Drin Diy 1A” (4.57) 
Ecs E¢sDrin Ecs4A cs 4 


Equation (4.57) shows that multiplying the parameter D~,, /€,, from Table 4.1 


by the desired D‘,, /A* and dividing the result by 4, the parameter D%,, /€,, for 


substituting D*,,/€ ., in Equation (4.56) is obtained. 


The same values of D%,, /A* correspond to the same and already derived values 


Of Peoy - 


D. PERFORMANCE ANALYSIS 
This Section revisits all performance analyses made in Chapter II by replacing 
the former pure AWGN scenario by a Nakagami fading channel. For this purpose, it 
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utilizes all results obtained in the two latter Sections concerning the estimation of P, for 
uncoded and coded M-QAM modulations as a function of the average signal-to-noise 
ratio per bit, Y,. The same three convolutional codes and four M-QAM modulators are 
presented. Each plot of P» versus Y,, is associated with a certain configuration and is 
shown for four values of m: 0.5, 1, 2, and 3. Recall that m = 0.5 represents the minimum 
value that m can assume and indicates a very severe fading condition, while m = 1 means 
a Rayleigh channel with no LOS, and m = 2 resembles a Ricean channel in which a LOS 
exists. Finally, m = 3 demonstrates a more favorable fading condition and can be 


considered a relative high value of m. 


1. 16-QAM with Dual-4, IEEE 802.11a and CHL 


Figures 4.4, 4.5, 4.6 and 4.7 show the performance of the 16-QAM in Nakagami 
fading channels for m = 3, 2, 1 and 0.5, respectively. In these figures, it is remarkable 
how large the coding gain becomes as the fading conditions worsen, reaching a value of 
more than 20 dB for m = 1 at P;, = 10°. It is interesting to notice how close to each other 
the uncoded curves become as m decreases. When m = 0.5, the uncoded 16-QAM, 8- 
QAM and QPSK curves unite themselves into a single curve. These three uncoded curves 
provide a means for comparing the performances of sending coded versus uncoded data 
through a channel that is operating at the bandwidth limit. Recall that the uncoded curve 
must be the QPSK when comparing the use of IEEE 802.11a and dual4 codes due to 
R. = 1/2 in both, while it changes to uncoded 8-QAM when considering the use of (4, 3) 
CHL. Although attractive in the pure AWGN scenario, the option of a smaller uncoded 
modulation becomes worthless in the Nakagami fading channel. Even considering a large 
m environment, as shown in Figure 4.4, notice that the uncoded QPSK curve does not 


interlace the coded curves as it did in the AWGN. 


Notice that the legends in Figures 4.4, 4.5, 4.6 and 4.7 display a percentage at the 
end of the CHL labels. This number indicates that an improvement in the upper bound 


was used for the correspondent curve. It represents the percentage of the total possible 


86 


symbol combinations that are covered by this improved upper bound, Po). Note, 
however, that if this percentage equals 100%, no improvements were applied. In order to 
avoid overloading the graphs with many curves that would lead to the same conclusions, 
the improvement technique was applied only for the best of the three codes, CHL, and 


limited to two steps. 


An interesting point to highlight is the maintenance of a constant step-gain for the 
improvements as the fading condition changes. As can be seen in Figures 4.4, 4.5, 4.6 
and 4.7, the average gains were about | dB and 2.5 dB for the two improvements and 


they remained constant in all four fading scenarios. 
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Figure 4.4. Probability of Bit Error for 16-QAM with SDD in Nakagami Fading 
(m = 3) Channel Using IEEE 802.11a, Duak4 and CHL Convolutional Codes. 
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Figure 4.5. Probability of Bit Error for 16-QAM with SDD in Nakagami Fading 
(m = 2) Channel Using IEEE 802.11a, Duak4 and CHL Convolutional Codes. 
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Figure 4.6. Probability of Bit Error for 16-QAM with SDD in Nakagami Fading 
(m = 1) Channel Using IEEE 802.1 1a, Dualk4 and CHL Convolutional Codes. 
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Figure 4.7. Probability of Bit Error for 16-QAM with SDD in Nakagami Fading 
(m = 0.5) Channel Using IEEE 802.11a, Dual 4 and CHL Convolutional Codes. 


In terms of relative performance among codes, the dual4 always slightly exceeds 
the IEEE 802.1 1a for all m’s. On the other hand, when comparing CHL with dual4 and 
IEEE 802.1 1a, it demonstrates that it is only better for m = 3. It is practically equal for 
m = 2 and loses approximately 2 dB when m is less or equal to 1. Therefore, the CHL 
seems to be more sensitive to fading than the other two. The reader should recall, 
however, that the CHL code presents the great advantage of compromising only 25% of 
the bit rate, whereas the use of either IEEE 802.11a or dua-k imposes 50% to this loss. 
Considering the pros and cons, the CHL code could also be considered the best choice 


among all three for 16-QAM in a Nakagami fading channel. 


2. 64-QAM with Dual-6 and IEEE 802.11a 


Figures 4.8, 4.9, 4.10 and 4.11 illustrate the performance of the 64-QAM 
associated with both dua-6 and IEEE 802.11a codes operating in four Nakagami fading 
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channels with m= 3, 2, 1 and 0.5, respectively. Notice again how large the coding gain 
becomes as the fading conditions worsen even using modulations with great M. Figures 
4.10 and 4.11 show coding gains greater than 20 dB at P; = 10° for both codes 


considered. 


The dual 6 code also proved to be better than IEEE 802.11a in Nakagami 
channels, ensuring its superiority already presented in the pure AWGN scenario. For this 
reason, the dual 6 code was plotted together with two improvements in the upper bound, 
which provided a 5 and 7.5 dB of gain for the channel analysis. Once more, these gains 
kept constant for the different m's, as can be seen in each one of the corresponding 
figures. Notice also that the step gains obtained for Nakagami channels due to 
improvements are the same achieved for AWGN, which can confirm the efficiency of 
this technique even for severe fading environments operating with large M modulation 


schemes. 
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Figure 4.8. Probability of Bit Error for 64-QAM with SDD in Nakagami Fading 
(m = 3) Channel Using IEEE 802.11a and Dual 6 Convolutional Codes. 
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Figure 4.9. Probability of Bit Error for 64-QAM with SDD in Nakagami Fading 
(m = 2) Channel Using IEEE 802.11a and Duat6 Convolutional Codes. 
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Figure 4.10. Probability of Bit Error for 64-QAM with SDD in Nakagami Fading 
(m = 1) Channel Using IEEE 802.11a and Dual-6 Convolutional Codes. 
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Figure 4.11. Probability of Bit Error for 64-QAM with SDD in Nakagami Fading 
(m = 0.5) Channel Using IEEE 802.11a and Dual 6 Convolutional Codes. 


3. BPSK and QPSK with Dual-k and IEEE 802.11a 


Figures 4.12, 4.13, 4.14 and 4.15 display the performance of the BPSK and QPSK 
when associated with both duaFk and IEEE 802.11a codes operating in four Nakagami 
fading channels with m = 3, 2, 1 and 0.5 respectively. Recall that the uncoded BPSK and 


QPSK, as demonstrated in Section IV.B, present the same performance. The differences 


in P;, come from the coding process that will provide different B (d)'s for each of these 
modulations. 


The IEEE 802.11a code proved to be better than dua-k for both BPSK and QPSK 
also in Nakagami channels, ensuring its superiority already established in the pure 
AWGN scenario. No improvements are possible to BPSK, as already seen in Chapter III, 
and one improvement is shown in all four figures for QPSK operating with IEEE 802.1la 
code. The same small improvement gain, approximately 0.7 dB, achieved in pure AWGN 


was again confirmed in the Nakagami channel actuating with QPSK. 
Q2 


Notice how great the IEEE 802.11a code performs with BPSK under severe 
fading conditions. In Figure 4.14, for m = 1, BPSK performance is even better than the 
one obtained for uncoded BPSK in pure AWGN. Also in Figure 4.14, IEEE 802.1 1a with 
QPSK and no improvement achieves almost the same performance as in pure AWGN for 
P;, = 10°. In other words, the IEEE 802.11a code is capable of eliminating a severe 
fading effect from the channel in both BPSK and QPSK at the expense of losing half 
bandwidth. 


Another interesting behavior to highlight is the close performance maintained 
between BPSK and QPSK through all different fading conditions when operating with 


duakk code. As the reader recalls, the duakk code keeps d,,, =4 for every M, which 


free 
greatly influences the close behavior previously cited. On the other hand, the IEEE 
802.1 1a code varies its djce according to Table 2.1. Since its dfee becomes smaller as M 


increases, the performance also decreases. 
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Figure 4.12. Probability of Bit Error for BPSK and QPSK with SDD in Nakagami 
Fading (m = 3) Channel Using IEEE 802.11a and Dual-k Convolutional Codes. 
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Figure 4.13. Probability of Bit Error for BPSK and QPSK with SDD in Nakagami 
Fading (m = 2) Channel Using IEEE 802.11a and Dual-k Convolutional Codes. 
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Figure 4.14. Probability of Bit Error for BPSK and QPSK with SDD in Nakagami 
Fading (m = 1) Channel Using IEEE 802.11a and Duak-k Convolutional Codes. 
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Figure 4.15. Probability of Bit Error for BPSK and QPSK with SDD in Nakagami 
Fading (m = 0.5) Channel Using IEEE 802.11a and Duakk Convolutional Codes. 


E. SUMMARY 
Chapter IV expanded the results originally obtained from the AWGN scenario to 
Nakagami fading channels. A second random variable expressing the fading effects now 


modeled a new upper bound of P,(d) for any M-ary coherent demodulator followed by 


SDD. The use of this upper bound in the special case of M-QAM modulations led to 
various performance analyses involving 16-QAM, 64-QAM, BPSK and QPSK under four 
fading figures m: 3, 2, 1 and 0.5. The reader had the opportunity to notice how the upper 
bound improvement technique is still valid and how convolutional codes can largely 
enhance performances under fading conditions. With Chapter IV, the last goals initially 
proposed for this thesis were achieved, leaving to the next chapter the task of 


summarizing all results. 
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Vv. CONCLUSIONS AND FUTURE WORK 


Many design tools and two analytical models were developed in this thesis. They 
have a general scope of applicability much larger than the special case of communication 
systems that make use of M-QAM modulations with Viterbi Soft-Decision decoding. The 
main conclusions of each of those tools and models will now be summarized per chapter 


together with a suggestion for future work at the end. 


A. CONCLUSIONS 

In Chapter II, a numerical method for determining the B (d) spectrum of any 
convolutional code was provided. This method turns into a general tool that can be 
applied to every convolutional code. As a specific application, it was shown how to 
adapt B (d) spectra obtained for binary convolutional codes to be used in many M-QAM 
modulation. The special case of the IEEE 802.1 1a code was analyzed and a table with its 
various B (d) spectra was presented. It is interesting to notice how its djree associated with 
each M-QAM modulation decays as M increases. For instance, the use of IEEE 802.11la 


convolutional code with BPSK leads to a d,,, =10, whereas this value drops to 3 when 


free 


the same code operates with 64-QAM. 


In Chapter III, a model for studying just the case AWGN channels was created. A 
simple formula for the upper bound of P,(d) for any M-ary coherent demodulator 
followed by a SDD was analytically derived. This formula is valid for M-QAM and also 
for the cases of orthogonal modulations such as M-FSK. A technique that improves the 
upper bound of P, (d) trying to tighten it to realistic data was shown. The possibilities of 
application and the efficiency of this technique were demonstrated depending on how 
large M is. As an example of extreme cases, no improvements are admitted to BPSK, 
while a realistic gain of up to 7.5 dB was achieved for 64-QAM. The performance 
analyses of 16-QAM, 64-QAM, BPSK and QPSK, associated with up to three different 
convolutional codes were shown. Those analyses demonstrated that the results for BPSK 


differ from QPSK when considering the different B (d) spectra they present for the same 
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code. The IEEE 802.11a code demonstrates itself better than dua-k code only when 
associated with BPSK and QPSK. When using 16-QAM and 64-QAM, the duakk code 
shows more efficiency. The CHL code demonstrates excellent performance with 16- 
QAM especially considering that it consumes only half the bandwidth consumed by dual 


k and IEEE 802.1 1a codes. Its coding gain reached 5.5 dB at P, = 10° when considering 


the use of improvements on the upper bound. 


In Chapter IV, the model for AWGN channels from Chapter III was expanded for 
the Nakagami fading channels, originating the second model of this thesis. A formula for 
the upper bound of P, (d) for any M-ary coherent demodulator followed by a SDD was 
analytically derived for Nakagami fading channels, showing itself valid for both M-QAM 
and for the cases of orthogonal modulations such as M-FSK. The performance analyses 
of 16-QAM, 64-QAM, BPSK and QPSK were repeated for the new model utilizing four 
fading figures m: 3, 2, 1 and 0.5. Considering all different configurations of modulations 
and codes, the same relative behavior among performances for the AWGN scenario was 
confirmed in Nakagami fading channels. As a general characteristic, convolutional codes 
associated to M-QAM with SDD demonstrated large coding gains, greater than 20 dB at 
P, =10°, under severe fading conditions. Also, this gain increases as the fading becomes 
worse. Concerning the upper bound improvements, they generate relative gains that 
maintain themselves constant as the fading figure of the channel varies. Additionally, as 
could be expected, those gains are equal to those provided in the pure AWGN scenario. 
The uncoded channels of the four studied modulations were also depicted and clearly 
demonstrated a trend for the same performances, no matter which M-QAM is being used, 


when operating in severe fading conditions (m <1.5 ). 


In the Appendix, a general formula to determine B (d) spectra of duakk codes is 


analytically derived. 
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B. FUTURE WORK 
Since the topics in this thesis have general application, many previous studies can 
now be revisited using the developed models in order to extend their results for the cases 


of M-QAM with SDD such as the performance analysis of OFDM in Reference [8]. 
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APPENDIX. DERIVATION OF £6(d) FOR DUAL-K 
CONVOLUTIONAL CODES 


The dua-k codes compound a class of nonbinary convolutional codes that are 
easily decoded by means of the Viterbi algorithm using either soft-decision or hard- 
decision decoding. It consists of two ( K = 2) k-bit shift-register stages and n = 2k 
function generators. Its output is two k-bit symbols for an input of one k-bit symbols 
(R. =1/2). The general form for the transfer function of a rate 1/2 duakk code is 
expressed as 
(*-)D' LN 


T(D,N, ) =——_———__.., 
1-NJ[2D+(2*-3)D?| 


where D and N are defined in Section I.B.4 and the exponent on J is equal to the number 


of branches in a given path [3]. 
The expression for T (D,N) can be obtained by making J = 1 in the formula for 
T(D.N,J), 


(2' -1)D*N 


ne) TN 2D +(28=3)D*] 


(A.1) 


For the first step, T(D,N ) will be expanded into an infinite series in powers of D 


and NV. Assume 
A=2? -1 = (A-2)=27=3, (A.2) 


and rewrite Equation (A.1) as 


AD*N 
TON)=—_—_—_\—. 
roe 1-N|2D+(A-2) D? | ve 
Let 
B= A-2, (A.4) 


so, Equation (A.3) can be developed as 


: AD*N Seer 2 
T(D,N) END = AD NUN 2D BD ]) 


co J co Jj 
= AD'NY"(2+BD) D'N! = AY\(2+BD) DN" 


j=0 j=0 


as j is i 
= B/E Jom ae pews 
j=0\ i= 


oo j y 
T(D,N)= ab E[/]p2™o jw 


j=0\ i=0\ 2 
Assuming that 
x= jt4+i => t=x-(jr+4), 
Equation (A.5) becomes 
oo ( 2j+4 . 
T(D,N)=A y BEG) I-44) Nit 
7a0| pea ES) 


Considering 


. = J (x-(j+4)) 4(2 j-x+4) 
c(Uj,x) = B Devt 
oe ae 


applied to Equation (A.7), the following can be obtained, 
co 2j+4 


T(D,N)=A) Y c(j.x)D*N“*? 


j=0x=j+4 


4 6 8 
T(D,N)= 4) nYco.no + vy c(1,x) D* + NY c(2,x)D* rm 
x=4 x=5 =6 


T(D,N)=A | Nc(0,4) D* 
+ N’ c(1,5)D° + N* c(1,6) D* 
+ N*c(2,6)D® + N*c(2,7)D’ + N*c(2,8)D* +... | 
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(A.5) 


(A.6) 


(A.7) 


(A.8) 


T(D,N) = AN c(0,4) D* 
+N°c(1,5) D” 
+ (N°c(1,6) + N’c(2,6))D° 
+ (N°c(2,7) + N*e(3,7)) D’ 


+ (N°c(2,8) + N*c(3,8) + N°c(4,8) ) D® 
+ ..] 
which yields 
°° d-3 
T(D,N)=A E| c(n-1, nD (A.9) 
d=4| n=f (d) 
d-2 
as (even d) 
where f(d) = a 
a (odd d). 


The second step is to find the partial derivative of Equation (A.9) with respect to 


N, and make N equal to | afterwards in order to obtain another series in powers of D 


only. 
oT(D,N lt ae 2 
ZBL ey nc(n-l, dyN°? | D4, 
oN d=4| n=f (d) 
hence 
dT(D,N = c % 
OT(DN) =p A n c(n-l, d) |D* =) B(d)D*. (A.10) 
ON |v deal nefa) d=4 





As can be seen, the expression for B(d) turns out to be the equation that defines 


the coefficients of this last series as a function of d. Hence, applying Equations (A.2), 


(A.4)and (A.8) to Equation (A.10) the following is obtained: 


B(d) = = OF = 1) y n lie = ie a B\o ee? qirn)-d +4) 
n=f (d) n— I)— 
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Finally, it can be concluded: 


d- 


—1 
B(d) = (2* -1) i aM : ,} = 3-9 (2n-d+2) 
n=f @) ~~ Pb 


w 


d=2 
2 
ged (odd d and d 25). 


(even d andd = 4) 
where f (d) = 


Notice that the minimum free distance djree = 4 symbols (4k bits). 
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